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"Would you tell me, please, which way I ought go from here?" 
"That depends a good deal on where you want to get to", said the Cat. 

"I don't much care where -" said Alice. 
"Then it doesn't matter which way you go, " said the Cat. 
"- so long as I get somewhere, " Alice added as an explanation. 
Oh, you 're sure to do that, " said the Cat, "if you only walk long enough. " 



from "Alice's Adventures in Wonderland" by Lewis Carroll 
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Introduction 



One of the main ideas of modern cosmology is that there was an epoch in the history 
of the Universe when vacuum energy dominated other forms of energy density such as 
radiation or matter. During this vacuum-dominated (or potential-dominated) era, the 
scale factor grew exponentially (or almost exponentially) in time and, for this reason, it is 
known as inflation. During inflation, a small and smooth spatial region of size less than 
the Hubble radius at that time could have grown so large to easily encompass the comov- 
ing volume of the entire Universe presently observable. If the early Universe underwent 
a period of so rapid expansion, one can then understand why the observed Universe is 
homogeneous and isotropic to such a high accuracy. All these virtues of inflation were 
already noted when it was first proposed in a seminal paper by Guth in 1981 [T]. Other, 
and more dramatic consequences of the inflationary paradigm were discovered soon after. 
Starting with a Universe which is absolutely homogeneous and isotropic at the classical 
level, the inflationary expansion will "freeze in" the vacuum fluctuations of the "inflaton" 
(the field responsible for the rapid expansion) which then becomes a classical quantity. 
On each comoving scale, this happens soon after the fluctuation exits the Hubble horizon. 
A primordial energy density, which survives after the end of inflation, is associated with 
these vacuum fluctuations and may be responsible for both the anisotropies in the Cos- 
mic Microwave Background (CMB) and for the large-scale structures of the Universe, like 
galaxies, clusters of galaxies, and dark matter. Inflation also generates primordial gravi- 
tational waves, which may contribute to the low multipoles of the CMB anisotropy. The 
amazing prediction of inflation is therefore that all of the structures we see in the present 
Universe are the result of quantum-mechanical fluctuations during the inflationary epoch. 
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Units and Conventions 



The so-called natural units, namely h—c—1, will be employed, so that 

[Energy] = [Mass] = [Length]"^ = [Time]~^ . 

nT'Pi — ^/VG is the Planck mass, whereas Mpi = l/y/SnG stands for the reduced Planck 
mass, and £p\ — y/Sn Mp^ is the Planck length. The present Hubble expansion rate 
is customarily parameterised as Hq = 100 /ikmsec"^ Mpc""*^ where h parameterises the 
uncertainty in the present value of the Hubble parameter. 

We shall also adopt the Einstein sum rule for repeated indices and the following 
conventions: The metric signature is (+,—,—,—), Greek letters denote space-time indices 
which run from to 3 and Latin letters denote spatial indices which run from 1 to 3. 

Some abbreviations used in this thesis are the following: 

CMB = Cosmic Microwave Background 

EW = elcctroweak 

BBN = Big Bang Nucleosintesis 

GR = General Relativity 

COBE = COsmic Background Explorer satellite 

DMR = Differential Microwave Radiometer 

WMAP = Wilkinson Microwave Anisotropy Probe 

ACBAR = Arcminute Cosmology Bolometer Array Receiver 

VSA = Very Small Array 

CBI = Cosmic Background Imager 

A CDM = ACold Dark Matter 

FRW= Fricdmann-Robertson- Walker 

WKB = Wcntzel-Kramers-Brillouin 

MCE = Method of Comparison Equation 

PS = Power Spectrum 

HFF = Horizon Flow Functions 

SR = Slow Roll 

HSR = Hubble Slow Roll 

PSR = Hubble Slow Roll 

GFM = Green Function Method 
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Units and Conventions 



Chapter 1 

What's the matter? 



What is the origin of the inhomogeneities we observe in the sky? What is their typical 
wavelength? How come that we are able to observe these inhomogeneities? The term 
inhomogeneity is rather generic and indicates fluctuations both in the background geom- 
etry and in the energy content of the Universe. These inhomogeneities can be represented 
by plane waves characterized by a comoving wave- number k. Since the evolution of the 
Universe is characterized by a scale factor a{t) which is a function of the cosmic time t, the 
physical wave- number is given by /Cph = k/a. The comoving wave-number is a constant 
and does not feel the expansion of the Universe, while the physical momentum is different 
at different epochs. Conversely the value of a given physical frequency is fully specified 
only by stating the time at which the physical frequency is "measured" . 

For example, if a given fluctuation has a momentum comparable with the present value 
of the Hubble parameter, A;ph,o — -ffo — 5.6x 10~^^ Mpi, we will have that A;ph,o — 2.3x10^^^ 
Hz. Fluctuations with momentum smaller that Hq have a wave-length larger than the 
present value of the Hubble radius and are therefore impossible to detect directly since the 
distance between two of their maxima (or minima) is larger than our observable Universe. 

At the decoupling epoch, when the radiation became free to propagate, the tem- 
perature of the Universe was of the order of a fraction of eV and the Hubble rate was 
i^dcc — 6.7 X 10~^^Mpi. The corresponding decoupling frequency was then red-shifted to 
the present value fcph.dec — 3 x 10~^^Hz. Fluctuations of this typical frequency appear 
as temperature inhomogeneities of the CMB induced by the primordial fluctuations in 
the background geometry and in the matter density. The latter fluctuations have already 
collapsed due to gravitational instability and formed galaxies and clusters of galaxies. At 
the epoch of the formation of light nuclei, the Universe was hotter, with a temperature of 
a fraction of MeV and the corresponding value of the Hubble expansion rate was about 
8 X 10~^^ Mpi. When the electroweak (EW) phase transition took place, the temperature 
of the Universe was of the order of 100 GeV leading, approximately, to i^EW ~ 10~^^ Mpi. 
Comoving wave-numbers of the order of the Hubble rate at the EW phase transition or 
at the epoch of the Big Bang Nucleosintesis (BBN) correspond, today, to physical fre- 
quencies of the order of A;ph,Ew ~ 2 x 10~^ Hz or fcph.sBN — 9 x 10~^^ Hz. Finally, in the 
context of the inflationary paradigm the value of the Hubble rate during inflation was 
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Chapter 1. What's the matter? 



Hini ~ 10~^ Mpi. Under the assumption that the Universe was dominated by radiation 
right after inflation, the corresponding physical frequency today is of the order of 10^ Hz. 

The detectors developed to measure anisotropics in the CMB are sensitive to frequency 
scales only slightly higher than 10^^^ Hz. Therefore, they cannot be used to "observe" 
fluctuations originated at earlier epochs and whose physical frequencies presently lie in 
the range 10~^ Hz < fcph < 1 MHz. For frequencies in this interval, one can only hope 
to detect the stochastic background of gravitational radiation which is therefore the only 
means we have to obtain direct information about the EW phase transition, the BBN and 
inflation itself. 

The cosmological inhomogeneities are usually characterized by their correlation func- 
tions. The simplest correlation function encoding informations on the nature of the fluc- 
tuations is the two-point function which is computed by taking averages of the fluctuation 
amplitude at two spatially separated points but at the same time. Its Fourier transform 
is usually called the power spectrum. On considering the physical scales discussed above, 
one is led to conclude that the power spectrum of cosmological inhomogeneities is defined 
over a huge interval of frequencies. For instance, the tensor fluctuations of the geome- 
try, which only couple to the curvature and not to the matter sources, have a spectrum 
ranging, in some specific models, from 10^^^ Hz up to about 1 GHz. 

From the history of the Universe sketched so far, some questions may be naturally 
raised. For our purposes, we just quote two of them: 

• Given that the precise thermodynamical history of the Universe above 10 MeV is 
unknown, how it is possible to have a reliable framework describing the evolution 
of the cosmological fluctuations? 

• Does it make sense to apply General Relativity (GR) to curvatures which can reach 
up to almost the Planck scale? 

As recalled above, we do not have any direct evidence for the thermodynamical state of 
the Universe for temperatures above the MeV, although some indirect predictions can be 
made. For example, the four lighter isotopes (D, ^He, ^Hc and ^Li) should have been 
mainly produced at the BBN below a typical temperature of 0.8 MeV when neutrinos 
decouple from the plasma and the neutron abundance evolves via free neutron decay. The 
abundances calculated in the homogeneous and isotropic BBN model agree fairly well with 
the astronomical observations. This is the last indirect test of the thermodynamical state 
of the Universe. For temperatures larger than 10 MeV, collider physics and the success of 
the standard model of particle physics offer a precious handle up to the epoch of the EW 
phase transition occurring for temperatures of the order of 100 GeV. There is then hope 
that our knowledge of particle physics could help us to flU the gap between the epoch of 
nucleosynthesis and that of the EW phase transition. 

The evolution of inhomogeneities over cosmological distances may well be unaffected 
by our ignorance of the specific history of the expansion rate. The key concept is that of 
conservation laws and, in the framework of metric theories of gravity (and of GR in partic- 
ular) , one can show that there exist conservation laws for the evolution of the fluctuations 
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which become exact in the hmit of vanishing comoving frequency k ^ 0. This conclusion 
holds under the assumption that GR is a valid description throughout the whole evolution 
of the Universe and similar properties also hold, with some modifications, in scalar-tensor 
theories of gravity (of the Brans-Dicke type) like those suggested by the low-energy limit 
of superstring theory. The conservation laws of cosmological perturbations apply strictly 
in the case of the tensor modes of the geometry whereas, for scalar fluctuations which 
couple directly to the (scalar) matter sources, the precise form of the conservation laws 
depends upon the matter content of the early Universe. Experiments aimed at the direct 
detection of stochastic backgrounds of relic gravitons which regard the high-frequency 
behaviour of cosmological fluctuations should therefore yield valuable clues on the early 
evolution of the Hubble rate. 

1.1 A short (incomplete) introduction to a long his- 
tory 

The history of the studies about the CMB anisotropics and cosmological fluctuations is 
closely related with the history of the standard cosmological model (see, for instance, the 
textbooks by Weinberg [2] and Kolb and Turner [3]). 

During the development of the standard model of cosmology, one of the recurrent issues 
has been the explanation of the structures we observe in the present Universe, which goes 
under the name of the structure formation. The reviews of Kodama and Sasaki and 
of Mukhanov, Feldman and Brandenberger ^ contain some historical surveys that are a 
valuable introduction to this subject. 

A seminal contribution to the understanding of cosmological perturbations is the paper 
of Sachs and Wolfe [6J on the possible generation of anisotropics in the CMB, where the 
basic general relativistic effects were estimated assuming that the CMB itself was of 
cosmological origin. The authors argued that galaxy formation would imply AT/T ~ 
10"^. Subsequent calculations [3 El El [TOl [H] brought the estimate down to 10~^, closer 
to what is experimentally observed. A second foundamental paper is the contribution 
of Sakharov [12] where the concept of acoustic oscillations was introduced. This paper 
contained two ideas whose legacy, somehow unexpectedly, survived until the present. The 
first idea is summarized in the abstract 

A hypothesis of the creation of astronomical bodies as a result of gravitational instability 
of the expanding Universe is investigated. It is assumed that the initial inhomogeneities 
arise as a result o/ quantum fluctuations... 

The model indeed assumed that the initial state of the baryon-lepton fluid was cold. So 
the premise was wrong, but, still, it is amusing to notice that a variety of models today 
assume quantum fluctuations as initial conditions for the fluctuations of the geometry at 
a totally different curvature (or energy). 

To summarise, Refs. [6] and [12] highlight three important problems regarding the 
fluctuations of the gravitational inhomogeneities: 
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Figure 1.1: Some of the most recent CMB anisotropy data from pU]: WMAP (filled 
circles); VSA (shaded circles) [20]; CBI (squares) [SHIS] and ACBAR (triangles) [23]. 

• their initial conditions; 

• their evolution; 

• their imprint on the CMB temperature fluctuations. 

1.2 CMB experiments 

Although we shall just study a particular aspect of the theory of cosmological perturba- 
tions, it is appropriate to give a swift account of our recent experimental knowledge of 
CMB anisotropies. A new twist in the study of CMB anisotropies came from the obser- 
vations of the COsmic Background Explorer satellite (COBE) and, more precisely, from 
the Differential Microwave Radiometer (DMR). The DMR was able to probe the angular 
power spectrum up to £ ~ 25 (see, for instance, [131 E] and references therein). 
DMR was a differential instrument measuring temperature differences in the microwave 
sky. The angular separation d of the two horns of the antenna is related to the largest 
multipole according to the approximate relation -i? ^ vr/f. 

COBE explored angular separations "i? > 7° and, after the COBE mission, various 
experiments have probed smaller angular separation, i.e., larger multipoles. A finally 
convincing evidence of the existence and location of the first peak in came from 
Boomerang [TSl IS] , Maxima [T7] and Dasi [IH] . These last three experiments explored 
multipoles up to £ = 1000, determining the first acoustic oscillation for £ ~ 220. An- 
other important experiment was Archeops [19] providing interesting data for the region 
characterizing the first rise of the angular spectrum. 

The angular spectrum, as measured by different recent experiments is reported in 
Fig. II. 1[ The data from the Wilkinson Microwave Anisotropy Probe (WMAP) are dis- 

^We just recall here that i{l + 1) Ciji^-K) measures the degree of inhomogeneity in the temperature 
distribution per logarithmic interval of t. Consequently, a given multipole i can be related to a given 
spatial structure in the microwave sky: small I will correspond to smaller and high i to larger wave- 
numbers around a present physical frequency of the order of 10~^^ Hz. 
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played as filled circles in Fig. 11.11 and provide, among other important evidences, a 
precise determination of the position of the first peak (i.e., i = 220.1 ± 0.8 [21]) and a 
clear evidence of the second peak. Further, once all the four-years data are collected, the 
spectrum will be obtained with good resolution up to £ ~ 1000, i.e., around the third 
peak. In Fig. 11.11 we also display the results of the Arcminute Cosmology Bolometer 
Array Receiver (ACBAR) as well as the data points of the Very Small Array (VSA) and 
of the Cosmic Background Imager (CBI). 

In recent years, thanks to the combined observations of CMB anisotropies [2S], Large 
scale structure [26l [27], supernovae of type la, BBN [28], some kind of paradigm for 
the evolution of the late time (or even present) Universe has emerged. This is normally 
called ACold Dark Matter (A CDM) model or, sometimes, "concordance model". The 
term ACDM refers to the fact that, in this model, the dominant (present) component 
of the energy density of the Universe is given by a cosmological constant A and a fluid 
of cold dark matter particles interacting only gravitationally with the other (known) 
particle species such as baryons, leptons, photons. According to this paradigm, our un- 
derstanding of the Universe can be summarized in two sets of cosmological parameters, 
one referring to the homogeneous background and one to the inhomogeneities. As far as 
the homogeneous background is concerned, the most important parameter is the indeter- 
mination on the (present) Hubble expansion rate which is customarily parameterised as 
Ho = hx 100km/[sec x Mpc] {h ~ 0.72 in the concordance model). There are various other 
parameters, such as the dark and CDM energy densities in critical units, flcdm — 0.224 
and ~ 0.73, respectively; the radiation and baryon energy densities in critical units, 

~ 8 X 10~^ and — 0.046; the number of neutrino species and their masses; the pos- 
sible contribution of the spatial curvature [Qk = in the concordance model) the optical 
depth of the plasma at recombination (of the order of 0.16 in the concordance model). 

The second set of parameters refers to the inhomogeneities. In the ACDM model the 
initial conditions for the inhomogeneities belong to a specific class of scalar fluctuations 
of the geometry which are called adiabatic and are characterized by the Fourier transform 
of their two-point function. This is the power spectrum which is usually parametrized in 
terms of an amplitude and a spectral slope. Besides these two parameters, one usually 
considers the ratio between the amplitudes of the scalar and tensor modes of the geometry. 

The purpose of this thesis is to report the development of a method for the resolution 
of the differential equation that drives the modes of scalar and tensorial perturbations 
of the cosmological inhomogeneities. This work is organized as follows: in Chapter [2] 
we review the Friedmann-Robertson- Walker (FRW) space-time and remark some aspects 
about gauge choices. Chapter [3] deals with the process of generation of cosmological 
perturbations by quantum fluctuations of a scalar field and a second-order differential 
equation for the quantities of interest is obtained. Such a kind of equations will be treated 
with semiclassical approximations in Chapter HI In Chapter E] we give new general results 
for the power spectra of cosmological perturbations, while in Chapter [6] we use those 
results to find, in detail, explicit expressions for the power spectra in several inflationary 
models. The thesis ends with our Conclusions and several Appendicies. 
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Chapter 2 



FRW Universes and their 
inhomogeneities 



2.1 FRW Universes 

A Friedmann- Robert son- Walker (FRW) metric can be written in terms of the cosmic time 
t as 



(0) 



1 — Kr^ 



+ (d^^ + sin^ ^d*^^) 



(2.1) 



where k can be 1, and —1 for spherical, eucledian or hyperbolic spatial sections, respec- 
tively; a is the scale factor, while r, 6, ip are the spherical coordinates. We can also write 
the metric in terms of the conformal time 1], defined by dr] = dt/a{t), as 



(0) 



g^ydx^ dx" = a {r]) < dr/ - 



dr^ 



+ {de^ + sin^ edip^) 



(2.2) 



The evolution of the geometry is driven by the matter sources according to the Einstein 
equations 



V = -R^ ,y — —6^,y R — 8 71 G T^i, , 



(2.3) 



where G^i, is the Einstein tensor, i?^,^ and R are the Ricci tensor and Ricci scalar, re- 
spectively; T^^ is the energy-momentum tensor of the sources. Under the hypothesis of 
homogeneity and isotropy, we can always write the energy-momentum tensor in the form 
diag (p, p, p, p) where p is the energy density of the system and p its pressure and 



are functions of the time. Using Eq. (12. 3p we obtain the Friedmann equations which, in 
conformal time, read 



= — — K 



n' 



AttG 



(p + 3p) 



(2.4) 
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Chapter 2. FHW Universes and their inhomogeneities 



along with the continuity equation 

p' + 3n{p + p) = . (2.5) 

In the above formulae 7i = a' /a and a prime denotes a derivation with respect to the 
conformal time. In cosmic time, the same equations become 



H' = ^-^--, H=-A7tG {p + p) + - (2.6) 

and the continuity equation is given by 

p + 3if(p + p) = 0, (2.7) 

where H = a/a is the Hubble expansion rate and dot denotes a derivation with respect to 
t. The relation between the Hubble factors in the two time parametrizations isTi. = aH. 



2.2 Fluctuations of the geometry in flat FRW Uni- 
verses 

Let us now consider a conformally fiat (i.e., k = 0) background metric of FRW type, 

9tiu{v) = a'^{v)Vtiu , (2.8) 
where 7]^^, is the fiat Minkowski metric. Its first-order fluctuations can be written as 

S gf,y{r], x) = Ss Qf^uiv, ^) + gfiuiv, ^) + Qf^uiv, ^) , (2-9) 

where the subscripts define, respectively, the scalar, vector and tensor perturbations. 
The perturbed metric S g^i, has ten independent components whose explicit form can be 
parametrized as 

6 goo = 2 a V 

6 Qij = 2a^ 5ij - di dj E) - a^ hij - a^ {di F, + d, Ft) 

Sgoi = -a^diB + a^Si, (2.10) 

with the conditions 

diS' = diF' = , hi = di hi = . (2.11) 

The decomposition given in Eqs. 02.101] is that normally employed in the Bardeen formal- 
ism [29] (see Ref. [5]). 

The scalar fluctuations are parametrized by the scalar functions (p, ip, B and E. 
The vector fluctuations are described by two divergenceless (see Eqs. (12. lip ) vectors in 
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three (spatial) dimensions Fj and Si (i.e., by four independent degrees of freedom). Fi- 
nally the tensor modes are described by hij, leading to two independent components (see 
Eqs. (EH])). 

Under general infinitesimal coordinate transformations 

T] ^ f] = T] + eQ , a;* X* = X* + e* (2.12) 

the fluctuations of the geometry defined in Eqs. f l2.10p transform as 

S g^,u Scjf,^ = 6 g^^ -V^e^ -Vu^f, , (2.13) 

where is the covariant derivative with respect to the background goemetry and = 
0^(77) (eo, —€i). These infinitesimal coordinate transformations form a group and the func- 
tions eo and are the gauge parameters. The gauge-fixing procedure amounts to fixing 
the four gauge parameters. Further, the spatial components can be separated into 
divergenceless and divergencefull parts 

ei = d^e + Q, (2.14) 

where di = 0. The gauge transformations involving eo and e preserve the scalar na- 
ture of the fluctuations while those parametrized by Q preserve the vector nature of the 
fluctuation. 

The background covariant derivatives in Eq. (12.131) are given by V^^ = e^— ^^^FJ^^ 
and, in terms of the unperturbed connections (see Eqs. (1A.7I) ). the fluctuations for the 
scalar modes can be written in the tilded coordinate system as 

0^0 = 0-Heo-eo (2.15) 

i^^i^ = tlj + neo (2.16) 

B^B = B + eo-e (2.17) 

E^E = E-e. (2.18) 

Under a coordinate transformation preserving the vector nature of the fluctuation, x* — >■ 
= X* + C (with di C = 0), the rotational modes of the geometry transform as 

S,^S, = S, + C (2.19) 
F,^F, = F,-Q. (2.20) 

Finally, the tensor fluctuations in the parametrization of Eqs. (12.101) are automatically 
gauge invariant. 

The perturbed components of the energy-momentum tensor can be written as 

5T° = 5p, 6T^ = -6p6i, 6T^ = {p + p) v , (2.21) 

where v is the velocity field and 6ui = di v. Under the infinitesimal coordinate transfor- 
mations (12.121) . the fluctuations given in Eq. (I2.2ip transform as 

Sp^Sp = Sp-p'eo (2.22) 
6p—^6p = 6p — wp'eo (2.23) 
V ^v = v + e' . (2.24) 
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Using the covariant conservation equation fl2.5p for the background fluid density, the gauge 
transformation for the density perturbation 6 = 6 p/p follows from Eq. fl2.22l) and reads 



We now have two possible ways of proceeding: 

• fix, completely or partially, the coordinate system (i.e., the gauge); 

• use a gauge-invariant approach. 
We shall briefly review both. 

2.3 Gauge fixing and gauge- invariant variables 

As we recalled above, the tensor modes are already gauge-invariant, whereas the vector 
and scalar modes change for infinitesimal coordinate transformations. Since the gauge 
fixing of the vector fluctuations is analogous, but simpler, to the one of the scalar modes, 
vectors will be discussed before scalars. 

2.3.1 Vector modes 

We have two different gauge choices for the vector modes. The first is to set Si = so 
that, from Eq. fl2.19l) . we have 



Since the gauge function is determined up to arbitrary (space-dependent) constants Cj, 
the coordinate system is not completely fixed by this choice. The second choice is = 
and, from Eq. (12.201) . Q = Fi. The gauge freedom is therefore completely fixed in this 
manner. 

Instead of fixing a gauge, we could work directly with gauge-invariant quantities. 
We note that, from Eqs. (I2.19P and fl2.2UI) . the combination F/ -|- Si is invariant under 
infinitesimal coordinate transformations preserving the vector nature of the fluctuations. 
This is the vector counterpart of the Bardeen potentials [29] . 

2.3.2 Scalar modes 

More options are available in the literature for scalar modes and we briefly review the 
more common ones (see Ref. [5] for more details). 




(2.25) 




(2.26) 
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Conformally Newtonian (or Longitudinal) gauge 

In this gauge we eliminate all the off-diagonal entries E and B of the perturbed metric in 
Eq. fl2.13p . i.e. Ei^ = B\^ = 0. From Eqs. fl2.17p and (12.181) the gauge functions are then 
determined as e = and eo = E' — B and the only non- vanishing entries of the perturbed 
metric are and ip. 

OfF-diagonal (or Uniform Curvature) gauge 

This gauge is characterized by the conditions Eod = and ipod = 0. From Eqs. (12.161) 
and (12.181) we find eg = —ip/Ti. and e = E. The two non- vanishing metric fluctuations are 
then B and (p. Since ipod = 0, the fluctuation of the spatial curvature also vanishes. 

Synchronous gauge 

This coordinate system is defined by 0s = i?s = 0, which however does not fix the gauge 
completely. In fact, Eqs. (12.151) and (I2.17p imply that the gauge functions are determined 
by 

(eoa)' = a0, e' = {B + e^) , (2.27) 

leaving, after integration over rj, two undetermined space-dependent functions. This gauge 
freedom is closely related with unphysical gauge modes. 

Comoving Orthogonal gauge 

In this gauge the quantity vc + Bq is set to zero and the expression of the curvature 
fluctuations coincides with a relevant gauge-invariant expression whose evolution obeys a 
conservation law. 

Uniform Density gauge 

With this choice the total matter density is unperturbed, 6 = 0. From Eq. (12.220 
the gauge function is then fixed as eo = {6 p)/p'. The curvature fluctuations on constant 
density hypersurfaces also obey a simple conservation law. If the total energy-momentum 
tensor of the sources is represented by a single scalar field one can also define the Uniform 
Field Gauge, i.e. the gauge in which the scalar field fluctuation vanishes 

Gauge-Invariant approach 

The set of gauge-invariant fluctuations can be written as [221 15] 



$ = + - EJ + n{B- E') 
^ = ^Ij -n{B - E') , 



(2.28) 
(2.29) 
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for the metric perturbations and 

Sp^ = 6p + p'{B-E') (2.30) 

6pg = 6p + p {B - E') (2.31) 

= v' + d'E' , (2.32) 

for the fluid inhomogeneities. Using Eqs. fl2.15p - fl2.18p and f l2.22l) - fl2.24l) . it can be verified 
that the quantities defined in the above equations are gauge-invariant. Eq. (12.321) can also 
be written as Vg = v + E' since = d^Vg and f * = 9* v. Sometimes it is also convenient 
to introduce the divergence of the pecuhar velocity 9 = div\ whose associated gauge- 
invariant quantity is 6 = (9j V^. 



Chapter 3 

Evolution of metric fluctuations 



We have seen in the previous Chapter that the tensor modes of the geometry are gauge- 
invariant. The vector modes are not invariant under gauge transformation, however, their 
description is rather simple if the Universe expands and in the absence of vector sources. 
The scalar modes are the most difficult ones: they are not gauge-invariant and couple 
directly to the sources of the background geometry. The analysis will be given for the 
case of a spatially fiat background geometry. 



3.1 Evolution of the tensor modes 



For the tensor modes the perturbed Einstein equations imply that 5tRi = 0. Hence, 
according to Eq. (]A.23p . we have in Fourier space 



/if + 2H/if + = 



(3.1) 



which is satisfied by the two tensor polarizations. In fact, since dih] = hi = 0, the 
direction of propagation can be chosen along the third axis and, in this case, the two 
physical polarizations of the graviton will be 



hi = -k 



ha 



hi 



h\ = hi^ 



(3.2) 



where h^ and h^ obey the same evolution Eq. (13.11) and will be denoted by h. Eq. (13.10 
can be written in two slightly different, but mathematically equivalent forms 







(3.3) 
(3.4) 



where /i^ = ah^. For k"^ 3> \a" /a\ the solutions of Eq. (13. 4p are oscillatory, and the 
amplitudes of the tensor modes decrease as 1/a, if the background expands. For k"^ -C 



^We recall that 5t, 5^ and 5s denote the first-order fiuctuation with respect to tensor, vector and scalar 
modes respectively. 
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\a"/a\ the solution is non-oscillatory and is given by 



dri' 

Hkiv) - ^ka{r]) + Bka{r]) / , (3-5) 



where and are numerical coefficients. 



3.1.1 Hamiltonians for the tensor problem 



The variable fik = ahk introduced above plays a special role in the theory of primordial 
fluctuations, since it is the canonical variable that diagonalises the action for the tensor 
fluctuations. Such a quadratic action for the tensor modes of the geometry can be obtained 
by perturbing the Einstein-Hilbert action to second-order in the amplitude of the tensor 
fluctuations of the metric. 

To this end, it is convenient to start from a form of the gravitational action which 
excludes, from the beginning, the known total derivative. In the case of tensor modes this 
calculation is rather straightforward and the result can be written as 



1 



gw 



64: nG 



(3.6) 



where rjap is the Minkowski metric. After a redefinition of the tensor amplitude through 
the reduced Planck mass 0, 



h = ^ — or h = ^ — 

the action (13. 6p for a single tensor polarization becomes 



(3.7) 



and its canonical momentum is simply given hj U = a'^ h'. The classical Hamiltonian 
associated with Eq. (13. 8p is 



d^x 



n 

— + {dihf 



(3.9) 



which is time-dependent. It is therefore possible to perform a time-dependent canoni- 
cal transformation, leading to a different Hamiltonian that will be classically equivalent 
to (13.91) . Defining the rescaled field fi = ah, the action (13.81) becomes 



Q(2) 
gw 



(3.10) 



^In our notation, Mpj^ = VSttG 
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which yields the Hamiltonian 



(3.11) 



for fi and of its canonically conjugate momentum n = fi' — Ti. fi. A further canonical 
transformation can be performed starting from (13.111) which is defined by the generating 
functional 



^2^3(/^,7r,r7) 



ax ( /X vr — — jj, 



(3.12) 



where n = fi' is the new momentum conjugated to fi. The new Hamiltonian can now be 
obtained by taking the partial time derivative of (I3.12p and reads 



(3.13) 



The possibility of defining different Hamiltonians is connected with the possibility of 
dropping total (non covariant) derivatives from the action. To illustrate this point, let us 
consider the action given in Eq. (13. 8p for a single polarization which can be rewritten as 



gw — I d''xd?7a^ h'"^ — {dihy 



(3.14) 



Recalling that, from the natural definition of the canonical field fi = ah, we also have 
ah' = ^' — H fi, Eq. (I3.14p becomes 



1 



S'gw = - / d xdr] jj! +H /i — 27Y/i/i' — {difxY 



(3.15) 



In Eq. (I3.15p . the canonically conjugate momentum is obtained by functionally differen- 
tiating the associated Lagrangian density C {rj, x), where 



-S'gw = J dr]L{r]) , L{r]) = j d^xC{r], 



x] 



(3.16) 



with respect to ^' and the result is n = ^' — Ti fi. Hence, the Hamiltonian will be 



Hgwiv) = J d^x [nfi' - Cg„ {r],x)] 



(3.17) 



Consequently, from Eqs. (I3.15P and recalling the expression of vr, one exactly obtains the 
result given in Eq. (13. lip . We now note that one can write 
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which can be replaced in Eq. fl3.15l) and leads to 



^ j d^x [/i'' + {n^ + W) fi^ - {d.fifj . (3.19) 



9 

The latter will be classically equivalent to the action fl3.15p for /i and its conjugate momen- 
tum vf = fi'. Hence, the Hamiltonian for fi and vf can be easily obtained from Eq. (13.171) 
with 71 replaced by tt and will exactly be the one given in Eq. (13.131) . 

The same conclusion also applies to the scalar fluctuations. In the Lagrangian ap- 
proach, this aspect reflects in the possibility of defining classical actions that differ by the 
addition of a total time derivative. While at the classical level the equivalence among the 
different actions is complete, at the quantum level this equivalence is somehow broken 
since the minimization of an action may lead to computable differences with respect to 
another. 

3.2 Evolution of the vector modes 

This paragraph is inserted for completeness, but in the following we only focus on the 
vector and scalar perturbations. 

The evolution of the vector modes of the geometry can be obtained by perturbing 
the Einstein equations and the covariant conservation equation with respect to the vector 
modes of the geometry 

<5^^- = 8 7rG(5vT; (3.20) 
V^5,T'^^ = 0. (3.21) 

Using Eqs. (1A.18I) and (1A.19I) of Appendix |Al the explicit form of the above equations 
can be obtained as 

V' S, = IQnG {p + p) a^V^ (3.22) 
Sl + 2nSi = (3.23) 
Vl + {l-3w)nVi = , (3.24) 

where is the divergenceless part of the velocity field and, as discussed in Chapter [21 the 
gauge has been completely fixed by requiring Fj = 0, implying, according to Eqs. (12.191) 
and (12.201) . that Q = Fi. Note also that, according to our conventions for the velocity 
field (see Eq. (jSi])), = {p + p) V\ Eqs. ^^-^Mj follow, respectively, from the 
— « and i — j components of the Einstein equations and from the spatial component of 
the covariant conservation equation. 

Eqs. fl3.22p and (13.231) are not independent. In fact, it is easy to check, using Eq. (13.221) . 
that Si can be eliminated from Eq. (13.231) and Eq. (I3.24p is then obtained. By Fourier 
transforming Eqs. (I3.22p and (I3.23p . we find the solution 

S.{ri) = ^ (3.25) 

V. iv) = ^'^f V , , (3.26) 
16 TT G a* [p + p) 
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where Ci{k) is an integration constant. Two distinct situations are possible. If the 
Universe is expanding, Si is always decreasing. However, Vj may also increase. Consider, 
for example, the case of a single barotropic fluid p = w p. Since, according to Eqs. (12.41) 
and (12.51) we have p ~ a~^^'^'^^\ then 

Vi oc a^""-^ . (3.27) 

U w = 1, Vi increases like r] (since 0(77) ~ y^)' "^^ile Si{ri) ~ r/"^ decays for large 
conformal times and is negligible for our purposes of describing the CMB. This observation 
was indeed raised in the past, in connection with the idea that the early Universe could 
be dominated by vorticity. If the Universe contracts in the future, the evolution of the 
vector modes will be reversed and Si may increase again. 

Finally, the situation may change even more radically if the theory is not of Einstein- 
Hilbert type or if it is higher-dimensional. Note that, if the energy-momentum tensor 
of the fluid sources possesses a non-vanishing torque, rotational perturbations can be 
copiously produced. 



3.3 Evolution of the scalar modes 

In analogy with the case of the tensor modes, large-scale scalar fluctuations follow a set 
of conservation laws that hold approximately for kt] <^ 1 and become exact only in the 
limit fci] — > 0. 

The evolution equations for the fluctuations of the geometry and of the sources can be 
written without fixing a specific coordinate system. Then, by using the definitions of the 
Bardeen potentials the wanted gauge-invariant form of the equations can be obtained. 
Inserting the definitions (I2.28I) - (I2.32I) into Eqs. flA.Sip and (IA.52p . the gauge-invariant 
Hamiltonian and momentum constraint become 

(^' + 7^<^>) =47ra^G5pg (3.28) 
(^' + 7^$) = -4 7rGa2 (p + p) e , (3.29) 

where B = diV^ is the divergence of the gauge- invariant (scalar) peculiar velocity field. 
Inserting the expressions of the gauge-invariant fluctuations given in Eqs. (I2.28I) - (I2.32I) 
into Eqs. flA.54p we can derive 



'^" + n ($' + 2^') + [n^ + ^n') $ + iv^($-^) =47rGa^6pg (3.30) 

o 

V2($-^) = 127rG'a^ (p + p) CT . (3.31) 

Finally, the gauge-invariant form of the perturbed covariant conservation equations be- 
come 

5Pg-3 {p + p)^' + {p + p) e + 3n {6pg + 6pg) = (3.32) 

(p + p) e'+[(p' + p')+4H (p + p)] Q + V^6pg + {p + p) V'$ = , (3.33) 
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as follows from inserting Eqs. fl2.28p - (l2.32p into Eqs. (1A.57P and (1A.58P without imposing 
any gauge condition. Eqs. fl3.28p and fl3.29p and Eqs. fl3.30p and fl3.3ip have the same 
form of the evolution equations in the longitudinal gauge. We also define the quantity 

n^-{^-nv,) = -^- )^,_!^, ^ , (3.34) 

which is gauge-invariant and will be used extensively later. 



3.3.1 Physical interpretation of curvature perturbations 

In the comoving gauge, the three-velocity of the fluid vanishes, i.e. vq = 0. Since 
hypersurfaces of constant (conformal) time should be orthogonal to the four-velocity, we 
will also have that Bq = 0. In this gauge the curvature perturbation can be computed 
directly from the expressions of the perturbed Christoffel connections bearing in mind 
that we want to compute the fluctuations in the spatial curvature, 

5 = 4 nf^ + 4 Rf^ ^ 4 ^'^c , (3.35) 

Gj 

where the subscript C is to remind that the expression holds on comoving hypersurfaces. 
The curvature fluctuations in the comoving gauge can be connected to the fluctuations in 
a different gauge characterized by a different value of the time coordinate, i.e. rjc ^ f] = 
?7c + cq. Under this shift 

V'c^^ = V'c + ^eo (3.36) 
{vc + Bc)^v + B= {vc + Be) + eo . (3.37) 

Since in the comoving orthogonal gauge vc + Bq = 0, Eqs. f l3.36p and ( I3.37P imply, in the 
new coordinate system, that 

^^ = ^^n {v + B) = -nVg) , (3.38) 

where the second equality follows from the definitions of gauge-invariant fluctuations 
given in Eqs (^M>-(^M>- From Eq. (Km . one can conclude that 7^ in Eq. (Km 
corresponds to the curvature fluctuations of the spatial curvature on comoving orthogonal 
hypersurfaces. 

Other quantities, defined in specific gauges, turn out to have a gauge-invariant inter- 
pretation. Take, for instance, the curvature fluctuations on constant density hypersurfaces 
ipi). Under infinitesimal gauge transformations we have 

V'D ^ = + ^ eo 

S pb ^ S p = 6 pb - p eo , (3.39) 

where, on constant density hypersurfaces, 5 Pd = by definition of the uniform density 
gauge. Hence, from Eq. f l3.39p . we obtain 

^B = ^ + n-j = <i/ + n-^ , (3.40) 
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where the second equahty follows from Eqs. (I2.28l) - fl2.32l) . Hence, the (gauge- invariant) 
curvature fluctuations on constant density hypersurfaces can be defined as 

C = -(^ + ^^) , (3.41) 

which coincides with the curvature fluctuations in the uniform density gauge. 

The values of ( and TZ are equal up to terms proportional to the Laplacian of This 
can be shown by using the definitions of TZ and (, whose difference gives 

C-n = -n (v, + ^-^^ . (3.42) 

From the Hamiltonian and momentum constraints and from the conservation of the energy 
density of the background 

[v^-^-3n {n^ + -^')] , (3.43) 



p' AnGa'^ {p + p) 
we in fact obtain 

^-^= 12.g!mp + p) - 

The density fluctuation on comoving orthogonal hypersurfaces can be also deflned as [29| 
_ 6pc _ 6p + p' jv + B) _ 5p^-3n jp + p) 
^""^ P P P ' 

where the second equality follows from the flrst one by using the deflnitions of gauge- 
invariant fluctuations. Again, using the Hamiltonian and momentum constraints, one 
obtains 

= ■ ^ , V^v^ , (3.46) 

which means that (^ — 7^) oc e^. 

3.3.2 Curvature perturbations induced by scalar fields 

Some of the general considerations developed earlier in this section will now be specialized 
to the case of scalar fleld sources characterized by a potential V{(p). For a single scalar 
fleld, the background Einstein equations In a spatially flat FRW geometry read 

n'=^-^ (C + VaA (3.47) 



3 V 2 

[n^ -n') =A7iG^'^ (3.48) 
dV 

^" + 2n^' + — a^ = . (3.49) 
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The fluctuation of ip deflned in Appendix |X] under a coordinate transformation (12.121) 
changes as 

X^X = X-^'(^o- (3.50) 

Consequently, the associated gauge-invariant scalar field fiuctuation is given by 

X = x + ^' {B-E') . (3.51) 

The two Bardeen potentials and the gauge-invariant scalar field fiuctuation X define the 
coupled system of scalar fiuctuations of the geometry, 



dip 



+ H ($' + 2^') + + $ = 47rG 



'2^1 / V/ 2 V 

—ip ^ + If X — — a X 

0(f 



(3.52) 
(3.53) 
(3.54) 



where Eqs. (13.521) . (13.531) and (I3.54p are, respectively, the perturbed — 0, — z and 
i — j components of the Einstein equations. Although these equations are sufficient to 
determine the evolution of the system, it is also appropriate to recall the gauge-invariant 
form of the perturbed Klein-Gordon equation 

d'^V dV 

X" + 2nX' - V'^X + —ra^ X + 2^ — a^ - p ($' + 3^') = . (3.55) 

The above equation can be obtained from the perturbed Klein-Gordon equation without 
gauge fixing, reported in Eq. ( ]A.50I) of Appendix lAj by inserting the gauge-invariant fiuc- 
tuations of the metric given in Eqs. (12^281) - (l2l29D into Eq. flA.501) together with Eq. fl3l5B . 

If the perturbed energy- momentum tensor has vansihing anisotropic stress, the {i 7^ j) 
component of the perturbed Einstein equations leads to $ = The gauge-invariant 
curvature fiuctuations on comoving orthogonal hypersurfaces are, for a scalar field source. 



7^ 



n(n<^ + ^') 



Eq. (I3.56P and a linear combination of Eqs. (I3.52p - p.54p lead to 



7^' 



n 



AnGp'' 



(3.56) 



(3.57) 



which implies that TZ is constant for the modes k such that 77 -C 1 [29] . The power spec- 
trum of the scalar modes amplified during the infiationary phase is customarily expressed 
in terms of TZ, which is conserved on super-horizon scales. 

Curvature perturbations on comoving spatial hypersurfaces can also be simply related 
to curvature perturbations on the constant density hypersurfaces. 



C 



m -n 



3 p 



(3.58) 
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Taking Eqs. (13.341) and (13.581) . and using Eq. (13.521) we have 



V?' 3^9 



/2 



2 ' 



(3.59) 



which confirms that TZ and C, differ by the Laplacian of a Bardeen potential. Taking the 
time derivative of Eq. (I3.57P and using Eq. (13.561) and Eqs. (I3.52I) - (I3.54I) . we obtain the 
second-order equation 



7^" + 2 - 7^' - v'^n = o , 



where the function z{r]) is 



z{ri) 



a if 



(3.60) 



(3.61) 



Eq. (I3.60p is formally similar to Eq. (13. ip . i.e. the evolution equation of a single ten- 
sor polarization. We can eliminate the first time derivative appearing in Eq. (I3.60p by 
introducing the gauge-invariant variable 



-zTZ 



(3.62) 



which leads to 



g - V^g = 



After a Fourier transformation, we have 

Qk + 

i.e., the scalar analog of Eq. (13. 4p . 



qk = 



(3.63) 



(3.64) 



3.3.3 Horizon flow functions 

The driving fields of the scalar and tensor modes are usually referred to as pump fields. 
They appear explicitly as the term a" /a in Eq. (13.40 and z"/z in Eq. (I3.64p . For general 
potentials leading to a quasi-de Sitter expansion of the geometry, they are conventionally 
described in terms of the hierarchy of horizon flow functions (HFF; see Appendix |B] for a 
comparison with different conventions), defined according to 

^° ~ H{N) 

din I e,- 1 , , 

e.+i = J>0 (3.65) 
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where N is the number of e-folds, N = ln(a/aj) [where a, = a{rii)] after the arbitrary 
initial time 77,. A general result is that inflation takes place for ei < 1. 

To find the expressions of the pump fields in terms of the HFF it is useful to write the 
background equations as a first-order set of non-linear differential equations 



- - 



dip J 2mpy 



ip = —2 nipi 



dH 

dip 



(3.66) 
(3.67) 



Bearing in mind Eq. (13.651) . the explicit relation determining the form of the pump field 
for tensor modes follows from 



W + n' = 2a^ + H = (2 - ei 



(3.68) 



where in the last equality, we use ei = —H/H^. In the same manner, for the scalar case 
we obtain 



^ 2 ..2 / 3 1 1 2 1 



(3.69) 



We conclude this section with a useful expression that relates 77, a, H and the HFF. 
Starting from the definition of the conformal time dr] = dt/a{t), we can integrate by parts 
and find that 



aH (1 -ei) 



(3.70) 



3.3.4 Hamiltonians for the scalar problem 

Different Hamiltonians for the evolution of the scalar fluctuations can be deflned. 

We begin by expressing the action of the scalar fluctuations of the geometry in terms 
of the curvature fluctuations 



S. 



(1) 

seal 



n'^ - {d{Rf 



(3.71) 



Using the canonical momentum -k-ji = z^TZ', the Hamiltonian related to the above action 
becomes 



TT. 



and the Hamilton equations 



z^ v^7^ , 7^' 



7^7^ 



(3.72) 



(3.73) 
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Eqs. f l3.73p can then be combined in the single second order equation (13.601) . 

The canonically conjugate momentum tt-ji is related to the density perturbation on 
comoving hypersurfaces. In particular, for a single scalar field source, one finds [29] 



p^ 



(3.74) 



where the second equality can be obtained using p^p + Pip = V^'^/a^ and the fact that the 
effective "velocity" in the case of a scalar field is = X/ip' . Inserting now Eq. (13.531) into 
Eq. (I3.52p . Eq. (I3.74p can be expressed as 



0? p^ 



3 IH? 



(3.75) 



where the last equality follows from the first equation in (12.41) . From Eq. (13.751) . it also 
follows that 



-6a Tier 



(3.76) 



in which we used Eq. (13.571) to express TZ' . Hence, in this description, the canonical field is 
the curvature fiuctuation on comoving spatial hypersurfaces and the canonical momentum 
is the density fiuctuation on the same hypersurfaces. 

To bring the second-order action in the simple form of Eq. (13.711) . various (non- 
covariant) total derivatives have been dropped. Hence, there is always the freedom of 
redefining the canonical fields through time- dependent functions of the background ge- 
ometry. In particular, the action (I3.7ip can be rewritten in terms of the variable q defined 
in Eq. ([322]). Then [30l [31] 



5, 



(2) 
seal 



d^X 



/2 



2 — qq 

z 



{^^q)' + 



(3.77) 



whose related Hamiltonian and canonical momentum are, respectively 



n = q 



In Eq. (I3.77P another total derivative term can be dropped, leading to the action 



S. 



(3) 
seal 



2 /d'x 



/2 \2 I ^ 2 

q - {diq) + —Q 



and Hamiltonian 



(3.78) 



(3.79) 



(3.80) 



where vr = q' . As in the case of the Hamiltonians for the tensor modes, Eqs. (13.720 . (13.780 
and (I3.80p are all related by canonical transformations. 
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3.4 Quantum amplification of metric fiuctuations 

In the framework of quantum theory, the evolution of the metric fluctuations can be 
described either in the Schrodinger formalism or the Heisenberg picture. 

In the Schrodinger description, the quantum evolution can be pictured as the spreading 
of a quantum mechanical wave-functional. The initial wave-functional is constructed 
as the direct product of states minimising the uncertainty relations for the harmonic 
oscillator modes forming the quantum field. This initial vacuum state has zero momentum 
and each Fourier mode of the field will evolve into a (two mode) squeezed quantum state. 
The uncertainty relations are still minimised on such states, but in such a way that one 
of the two canonically conjugate operators will have a variance much larger than the 
quantum limit, while the other canonically conjugate operator will have a variance much 
smaller than the quantum uncertainty limit. As a background pumping electromagnetic 
field (a laser) is able to produce, under some circumstances, squeezed states of photons, 
the classical gravitational field (for instance the curvature) is able to produce squeezed 
states of gravitons (for the tensor modes) or of phonons (in the case of the scalar modes). 

In fully equivalent terms, the evolution of the fluctuations can be described in the 
Heisenberg representation. This is the description which will be adopted here. The 
expectation values of the Hamiltonians of the scalar and tensor modes of the geometry 
will then be minimized for r] —>■ —oo, which is physically identified with the beginning of 
inflation. 

3.4.1 Large-scale power spectra of tensor fluctuations 

Let us consider the Hamiltonian given in Eq. (13.131) and drop the tildes in the momenta. 



where fi = a Mpi/ (an identical expression holds for the other polarization). After 
imposing the commutation relations for the canonically conjugate quantum fields (in 
natural units with h = 1) 




(3.81) 



[fi{x,r]),Tr{y,r])] =i5^^\x - y), 



(3.82) 



the operator corresponding to the Hamiltonian fl3.8ip becomes 




(3.83) 



In Fourier space the quantum fields can be written as 




(3.84) 
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Inserting now Eqs. (13.841) into Eq. fl3.82p and demanding the validity of the latter implies 
the following canonical commutation relations for the Fourier components of the quantum 
operators 



(3.85) 



Inserting now Eq. (13.841) into Eq. (I3.83P we get the Fourier space representation of the 
quantum Hamiltonian|j 



(3.86) 



H{v) = 4 I d'k 



In the Heisenberg representation the field operators obey 



il,H 



and i if' 



Using now the mode expansion (I3.84p and the Hamiltonian ( 13T 
Fourier components of the operators is determined by 



whose general solution is 

where the mode function fk obeys 

n + 



«fc(^o) fkiv) + al^Vo) fkiv) 



(3.87) 

the evolution for the 

(3.88) 
(3.89) 

(3.90) 
(3.91) 



fk = 0, 



(3.92) 



and gk = /(. In Eqs. fl3.90p and (13.910 the operators a^(?7o) annihilate the vacuum, 



%(?7o)|0_oo) = for ?7o 



-oo . 



(3.93) 



As anticipated, this is the state which minimises the Hamiltonian (13.830 close to the onset 
of inflation. We remark that the time t]o is the same for all the modes k. The second 
remark is that the vacuum has zero total momentum during the entire evolution. This 



■^Note that, in order to derive the following equation, the relations fi^ - = fij^ and tt^ - = tt^ should be 
used. 
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means that particles can only be produced in pairs from the vacuum triggered by the 
pumping action of the gravitational field. Inserting Eqs. (l3.9Up and fl3.91l) into Eq. (13.841) . 
the following Fourier expansions can be simply obtained 



fi{x,r]) 
fc{x,ri) 



(2 7r)3/2 
1 

(2 7r)3/2 



d^A; 



i k-x 



(3.94) 
(3.95) 



since, upon integration by parts, one has a_j:e 



i k-x 



-i k-x 



and d ^ e 



t „—ik-x "t„jfc- 



The two point function for the tensor fluctuations can then be obtained as 



k 



0_ 



h{x,7]) h{y,r]) 



0. 



2/2 



d'k\Mv)\'e 



2 —i k-f 



(2 7r)3a2(r7) 

where r = {x — y). After angular integration, the previous expression becomes 

sin kr 



(3.96) 



0_ 



0- 



dlnkVh 



kr 



(3.97) 



where 



Vn = ^^-,\fk{v)f 



(3.98) 



is the tensor power spectrum. In the study of inflation, the tensor spectral index is usually 
defined as 



and its a-running is 



dlnPfe 
dln/c 



(dlnA;)2 ■ 



(3.99) 



(3.100) 



3.4.2 Large-scale power spectra of scalar fluctuations 

The calculation of the scalar power spectrum follows exactly the same steps discussed in 
the case of the tensor modes. 

In full analogy with the calculation of the tensor power spectrum, the Hamiltonian 
given in Eq. (13.801) can be chosen. Promoting the canonical mode q and its conjugate 
momentum to quantum operators which satisfy 



[q{ri,x),'kq{ri,y)]=i5''^\x-y) , 



(3.101) 
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one can also define a vacuum state |0_oo) for the scalar modes at rj ^ — oo. Recalling 
the relation between the canonical modes and the gauge-invariant curvature fluctuations, 
obtained in Eq. (13.621) . we can then write the scalar two-point function as 



0_ 



0_ 



(3.102) 



where f^j^{f]) are the mode functions pertaining to the scalar problem and obeying the 
equation 



s k 



. 



(3.103) 



The scalar power spectrum, the scalar spectral index and its running will be deflned, with 
a convenient redeflnition of the tensor quantities, in the Chapter O 
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Chapter 4 

WKB-type approximations 



In this chapter we give a review of some Wentzel-Kramers-Brillouin (WKB)-type ap- 
proaches used to solve second-order differential equations. We start in the first section by 
reviewing the standard WKB approximation and then proceed, in the second section, by 
presenting the problem that arises around the turning point (s) of the frequency and how 
it can be avoided in other WKB-type approximations. In particular, we give a first im- 
proved version of the WKB analysis in the third section and a generalized WKB method, 
also called Method of Comparison Equations, in the last section. 

4.1 Standard WKB approximation 

In this Section we begin by recalling the standard WKB method, as is introduced in 
almost all Quantum Mechanics textbooks (See, for example [H^). 

The WKB approximation is defined by first introducing a "small" parameter 6 > 
in order to perform the formal adiabatic (or semiclassical) expansion of the mode func- 
tions [33]. Let us briefly recall that such an expansion is called adiabatic when \/6 is 
the inverse of a (typically very large) time over which the system evolves slowly, and 
semiclassical when ~ h. Further, the limit 6^1 must be taken at the end of the 
computation. Consequently, second-order differential equations such as Eqs. 03.41) and 
f l3.64p [or equivalently Eqs. fl3.92p and (13.1031) ] are formally replaced by 



■ d^ 1 

S — + n^{k,r]) /i = 0. 



(4.1) 



We then denote the leading order term in such an expansion of /i by 



(4.2) 




which satisfies the following differential equation 



S— + Q'^{k,r])-SQn{k,r]) /i^KB = , 



(4.3) 
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where 



Qnik,r]) 



3 jn'y 

4 ^]2 



2Si 



(4.4) 



We now observe that /xwkb is the exact solution of Eq. (14.1 1) in the "adiabatic hmit" 
6^0, and is expected to be a good approximation to the exact /i (in the hmit 6 —>■ 1) if 



A 



fi2 



< 1 . 



(4.5) 



Sometimes, to respect the above condition, we need to perform a suitable change of 
variables. This will indeed be the case for the cosmological perturbation and in Chapter [5] 
we shall give the correct transformation to improve all WKB-type approximations. After 
that change of variables we shall use the following form for the differential equation 



dx^ 



+ u;^(x) 



X 







(4.6) 



which is itroduced here to keep track of that transformation. Finally, we recall that 
Eq. (14. 6p . but also Eq. (14. ip . in general suffer of some problems near the turning point of 
the frequency a;(x) which we now discuss. 



4.2 The "Turning Point" problem 

We describe here a kind of problem that can arise around the turning point of the fre- 
quency for the standard WKB approximation and how we can avoid it in other WKB-type 
approximations. First of all, a "Turning Point" is defined as that point where the fre- 
quency of a second-order differential equation like (14. ip or (14. 6p vanishes (see Fig. 14. ip . 
A solution of the type given in Eq. (14. 2p (or the analogue for Eq. (14. 6p ) diverges in this 
point (as well as the error given by the term in Eq. (14. 4p ) and it is well known that one 
therefore needs a different branch of solution there. 

Schematically, for the standard WKB approximation we have the situation shown in 
Fig. 14.21 where we see the solution for an ordinary second-order differential equation with 
only one turning point. The standard procedure to address this problem is to divide the 
domain of the independent variable into three regions: on the left and on the right of the 
turning point, where the standard WKB solution is valid, and the "turning point region", 
where we use a non-WKB solution. Usually, around the turning point, the true frequency 
is then replaced by a linear approximation and the approximate solution is then given by 
a linear combination of Airy functions. 

The next step is to connect the three branches to obtain a "global" solution on the 
whole range of the independent variable. We usually match the WKB solutions (i.e., the 
oscillatory and exponential branches) to that around the turning point using two asymp- 
totic expansions of the latter branch. The overlap regions, where this asymptotic match is 
performed, are not a 'priori defined but should be determined so as to minimise the error. 
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-1 -0.5 0.5 1 



Figure 4.1: Sketch of a typical frequency with one turning point (as is the case for cosmo- 
logical perturbations). The independent variable is irrelevant for the discussion developed 
here. 

which is usually a very difficult task, although it does not normally affect the accuracy 
of the solution to leading order in 5 (in Fig. 14.21 we use two grey bands to indicate their 
possible location). This subtlety has already been studied by the physicists and mathe- 
maticians who first used the WKB approximation and it is known that the method can 
give some problems, for example, if one needs a more accurate solution at next-to-leading 
order in the expansion in 5. In fact, if we used next-to-leading order expressions of the 
standard WKB expansion, we could find that the overlap regions, where the match is to 
be performed, become too narrow (or disappear at all!) and an accurate match of the 
different branches is no more possible [Ml ES] • This could even lead to next-to- leading 




-0.5 0.5 1 1.5 



Figure 4.2: Standard WKB approximation. 
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order results less accurate than the leading ones [SSj . 




-0.5 0.5 1 1.5 



Figure 4.3: Improved WKB approximation. 

The improved WKB approach, that we shall present in detail in the next section, will 
give us an approximation which is also valid at the turnig point. The solution will be 
divided into only two branches and the procedure of matching will be clearer and more 
accurate. Actually, the match will be performed exactly at the turning point, which is 
well located, and, without ambiguous matching regions, the results at next-to-leading 
order are always improved (see Fig. 14.31) . 

Finally, the Method of Comparison Equations, presented in the last part of this chap- 
ter, is a generalised WKB method that will not need any matching between branches 
of solutions. The solution will be valid on the whole range of the independent variable 
and could be used in every problems where we need a global solution of our differential 
equation (see Fig. 14.41) . 

So far we haven only considered one turning point, since this is the case for the 
equations governing cosmological perturbations we shall solve for later. We however note 
that, if the frequency has two, or more, turning points the solution of this problem could be 
enormously complicated (a case with two turning points will be described in Appendix [Gj. 

4.3 Improved WKB approximation 

We shall now describe a method to improve the WKB approximation for equations of 
the form in Eq. (14. 6 p to all orders in the adiabatic parameter 6 by following the scheme 
outlined in Ref. |37j (see also Ref. [38]). The basic idea is to use approximate expressions 
which are valid for all values of the coordinate x and then expand around them. Two 
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Figure 4.4: MCE approximation. 

different such expansions we shall consider, namely one based on the Green's function 
technique and the other on the usual adiabatic expansion [37j . 

Let us first illustrate the main idea of Ref. [37] for the particular case u"^ = C (x — xq)", 
with C a positive costant. All the solutions to Eq. (14. 6 p for x > xo can be written as 
linear combinations of the two functions 



U±[X) 



Lu{x) 



(4.7) 



where 



^(x) = / uj{y)dy , m 



XQ 



n + 2 



and Ji, are Bessel functions [o9]. Moreover, for a general frequency u, the expressions in 
Eq. diZD satisfy 



dx^ 



+ uj'^{x) — cr(x) 



u±=0 



(4.9) 



where the quantity (primes will denote derivatives with respect to the argument of the 
given function from here on) 



3 (uj 



i\2 



a[x] 



(4.10) 
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contains the term Q defined in Eq. ( 14.41) . whose divergent behavior at the turning point 
X = Xq we identified as the possible cause of failure of the standard WKB approach. For 
a general (finite) frequency, which can be expanded in powers of x — xq 0, 



uj'^{x) = C (a; — XqY 



9>1 



(4.11) 



one finds that the extra term in Eq. fl4.1Up precisely "removes" the divergence in Q at the 
turning point. In fact, the residue 



3C2 



. X 3 (n + 5) cl 

lim aix) = — ; r 

x^xo ^ ^ 2 (n + 4) (n + 6) n + 6 



(4.12) 



is finite The finiteness of a at the turning point is crucial in order to extend the 

WKB method to higher orders. For slowly varying frequencies, with |cr(a;)| small every- 
where, the functions (14.71) are thus expected to be good approximations to the solutions 
of Eq. (14. 6 p for the whole range of x, including the turning point. 

We can now introduce both the adiabatic expansion of the previous Section and a new 
formal expansion of the mode functions by replacing Eq. (14. 6 p by 



6^ ■ 
dx^ 



+ uj (x) — 6 a{x) 



X = -5ea{x)x 



(4.13) 



where S and e are "small" positive parameters. For clarity, we shall refer to expressions 
proportional to 5" as the n-th adiabatic order and to those proportional to as the n-th 
perturbative order. 

It is convenient to consider the solutions to Eq. (I4.13P to the left and right of the 
turning point Xq separately. We shall call: 



region I) where cj^ > (on the right of xq), with 

uji{x) = \/u'^{x) , ^i(a;) e 



(4.14a) 



Xo 



and the solutions to the corresponding homogeneous equation (14. 9 p are given in 
Eq. (S7D; 



region II) where u;^ < (on the left of xq), with 

LJu{x) = ^/-LJ^{x) , ^11 (x) 



xo 



t^ii(y) 



(4.14b) 



and the solutions to the corresponding homogeneous equation (14. 9p are obtained 
from those in Eq. (14. 7p by replacing the J^, with the I^, [39j . 



^In the following we shall consider the general case in Eq. (|4.1ip . although one might restrict to n = 1, 
since the frequencies for cosmological perturbations usually exhibit a linear turning point. 



4.3. Improved WKB approximation 



51 



Corresponding expressions are obtained for all relevant quantities, for example ai and an, 
and we shall omit the indices I and II whenever it is not ambiguous. 

Although it is possible to expand in both parameters, so as to obtain terms of order 
6^ e'^ (with p and q positive integers), in the following we shall just consider each expansion 
separately, that is we shall set 5 = 1 in the perturbative expansion (see Section 14.3.11) and 
e = 1 in the adiabatic expansion (see Section I4.3.2p . 



4.3.1 Perturbative expansion 

For 5 = 1, the limit e — > is exactly solvable and yields the solutions (14.71) . whereas the 
case of interest (14.61) is recovered in the limit e 1, which must always be taken at the 
end of the computation. The solutions to Eq. (14.131) can be further cast in integral form 

as 



xix)=u{x)-6 / G{x,y)a{y)xiy)dy 



(4.15) 



where u{x) = Aj^u+{x) + A„m_(x) is a linear combination of solutions (14.71) to the 
corresponding homogeneous equation (14.91) . and G{x, y) is the Green's function determined 

by 



dx^ 



+ uj'^{x) — cr(x) 



G{x,y) = e5{x-y) . 



(4.16) 



The solutions of the homogeneous equation (14. 9 p in the two regions are then given by 

ui{x) = 



UJiiX 



(4.17) 



Unix 



UJuix) 



l^uix)] + B_ [^u{x)]} = B+uu+{x) + B_ un^{x) 



and the Green's functions can be written as 

TT TT 

Gi{x,y) = --j=e{x-y)ui^{x) ui+{y) - -j=9{y - x) ui+{x) ui_{y) 



t TT % TT 

Gii{x,y) = -^6{x - y) uii^{x) uii+{y) + —9{y - x) un+{x) uu-iy) ■ 
With the ansatz 

Xi{x) = a+(x) Mi4.(x) + y4_ a_(x) Mi_(a;) 
Xii(x) = B+b+{x)uii+{x) + B_b-{x)uii-{x) , 



(4.18) 



(4.19) 
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from Eqs. fl4.17p and fl4.18p . the problem is now reduced to the determination of the 
x-dependent coefficients 



a+{x) = [l + e J+_(x, Xi)~\ + e —— J 



IX, Xi] 



(4.20) 



a_(x) = [l +e J^_(xo,x)] + £ -j— -'^J+(a;o, a;) 

b+{x) = [l - ie I^_{x,xo)] - ie — IZ^{x,xq) 

b-{x) = [l - ie I^_{xf,x)] -ie^lX+{,Xf,x) , 

in which, for the sake of brevity (and clarity), we have introduced a shorthand notation 
for the following integrals 

•'X2 



Jss{^i,X2) = -y= J cTi{y) a^{y) uis{y) uis{y) dy 



Vs Jxi 

[ ^^~t\ ^«'(^) ^^"^ [^i(^)] "^^m My)] dy 



Iss{xi,X2) = -j= / aii{y) bn,{y) uiis{y) uus{y) dy 



(4.21) 



"\/3 Ai 

^ . ^ eii(i/) 

V3 ^11 (y) 



&«,(y) /sm [^Il(l/)] ^Sm [^u{y)] dy , 



where the symbols w, s and s are + or — . 

Before tackling this problem, we shall work out relations which allow us to determine 
the constant coefficients A± and B± uniquely. 

4.3.2 Adiabatic expansion 

Let us now apply the usual adiabatic expansion with the assumption that the leading 
order be given by the functions (14.71) [rather than the more common expression (14. 2p ]. 
This leads one to replace 

e(x)^^ = 4^ rcu{y)dy, (4.22) 



V6 V6 
and consider the forms [?7] . 

U± = F±(x; 5) u±{x) + G±{x; 5) u'^{x) , (4.23a) 

where the x-dependent coefficients F and G must now be determined. The particular 
case of a linear turning point [30] (precisely the one which mostly concerns us here), can 
then be fully analyzed as follows. 
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Recursive (adiabatic) relations 

The x-dependent coefficients F and G can be expanded in powers of 6 as 

F(x;5) = ^5^>0)(x) 

i>o 

(4.23b) 

G{x;6) = ^6^^^j){x) . 

Upon substituting into Eq. (14.131) with e = 1, one therefore obtains that the coefficients 
and 7(j) are given recursively by the formulae 



1 r 

^^^■^^"^^ " 2Mx)J ^""^^^ [2 7|._i)(y) + 0(,-_i)(y)] 

+(^'iy) i(j-i)iy) + (l>{j^i)iy)} 



(4.23c) 



where it is understood that the integration must be performed from xq to x in region 1 
and from x to Xq in region II. Let us also remark that, for u;^ ~ (x — Xq) near the turning 
point, the above expressions are finite [lOj, whereas for more general cases one expects 
divergences as with the more standard WKB approach. 

Leading order 

The lowest order solutions (14. 7p are correctly recovered on setting 

0(0) = 1 and 7(0) = , (4.24) 

in the limit 5 = 1. The relevant leading order perturbations are therefore obtained on im- 
posing the initial conditions (15.131) and matching conditions between region I and region II 
at the turning point, which therefore yield for xi and xii the same linear combinations ui 
of Eq. (I5.25ap and Uu of Eq. fl5.25bl) previously obtained. 

Higher orders 

The condition (14.241) makes the formal expressions of the coefficients and 7(1) partic- 
ularly simple. 



= I b(i)(^) + ^y 



7(1) a; = , . / ^^dy 



(4.25) 
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and the first order solutions are then given by linear combinations of the two functions 
U±ix)= + M±(a;)+57(i)(a:)u'±(a;) . (4.26) 

4.4 Method of Comparison Equation 

In this Section we briefly review the Method of Comparison Equation (MCE). The name 
is due to Dingle and it was independently introduced in Refs. [HI US] and applied to wave 
mechanics by Berry and Mount in Ref. [13]. The standard WKB approximation and its 
improvement by Langer [37] are just particular cases of this method. In Section [4.4. II we 
shall also present the connection of the MCE with the Ermakov-Pinney equation. 

We start from the well-known second-order differential equation 04.61) and suppose 
that we know an exact solution to a similar second-order differential equation, 

"d^ 



da2 



U{cr) = , (4.27) 



where G is the "comparison function". One can then represent an exact solution of 
Eq. (14.61) in the form 

Xix) = (^] Uia) , (4.28) 



dx / 

provided the variables x and a are related by 

Eq. fl4.29p can be solved by using some iterative scheme, in general cases as [44J and 
Appendix [C] or Refs. [13, Il6] for specific problems. If we choose the comparison function 
sufficiently similar to a;, the second term in the right hand side of Eq. (14.291) will be 
negligible with respect to the first one, so that 

u\x) ^ (^^^ 02((t) . (4.30) 

On selecting a pair of values Xq and ctq such that ctq = cr(xo), the function a{x) can be 
imphcitly expressed as 

- e(x) = r v/±^dy ^ r v/±0%^dp , (4.31) 

•J xn J cm 



where the signs are chosen conveniently %. The result in Eq. (14.281) leads to a uniform 
approximation for xi^)-, valid in the whole range of the variable x, including "turning 
points" . The similarity between B and uj is clearly very important in order to implement 
this method. 



^Here ^(x) corresponds to the quantity Cii(a;) used in Section l¥751 when the signs are both chosen as 
minus. 
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4.4.1 MCE vs Ermakov&Pinney 

In Ref. |17] , we studied the relation between the MCE l^j and the Ermakov - Pinney 
equation [HI HQ] . Let us consider again the second-order differential equation (14. 6 p but 
with the small parameter e defined by [f| 

^ = -yix)x{x) . (4.32) 
Let us suppose that one knows the solution of another second-order differential equation 

^^eVM.). (4.33) 

In this case one can represent an exact solution of Eq. (I4.32p in the form 
where the relation between the variables x and a is given by the equation 



dx J \dx J dx^ \dx J 

On solving Eq. (14.351) for a{x) and substituting for it into Eq. (14.341) one finds the solution 
to Eq. (14.321) . Eq. (I4.35P can be solved by using some iterative scheme. The application 
of such a scheme is equivalent to the application of the MCE or, in other words, to 
the construction of the uniform asymptotic expansion for the solution of Eq. (I4.32p . The 
method of construction of the solution to Eq. (I4.32p by means of the solutions of Eqs. (14.330 
and (I4.35P is called the method of comparison equations and the function G^(o") is called 
the comparison function [32| . 

The iterative scheme of solution of Eq. (I4.35P depends essentially on the form of the 
comparison function S(ct). A reasonable approach consists in the elimination of a from 
Eq. (14.350 and its reduction to a form, where the only unknown function is (da/dx). The 
equation obtained for the case of the simple comparison function ©^(cr) = 1 coincides with 
the Ermakov-Pinney equation, this will be explicitly shown in the next section. For more 
complicated forms of the comparison function 0^((t) useful for the description of various 
physical problems the corresponding equation will acquire a more intricate form and we 
shall call it a generalised Ermakov-Pinney equation. The role of an unknown function in 
this equation is played by the function 



•^With respect to Eq. ()4.6p . we have absorbed a minus sign in e, as this change of sign is irrelevant for 
the present topic and for its apphcation in Appendix [Cl 
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while the equation itself has the form 

A-(,.^)=.l(,.^.^.^) , (4.37) 

\ ax J \ ax dx^ ax-^ J 

where the explicit forms for the functions K and L will be defined in the next sections. 
We shall represent the function y{x) as a series 



y = J2yne'', (4.38) 



n=0 

and the general recurrence relation connecting the different coefficient functions i/n can 
be represented as 

= . 4.39 

de- de^-^J 

To implement this formula one may use the combinatorial relation 

d'^iuv) , („) 



dx^' 



(u + ^;)^"^ (4.40) 



where on the right-hand side one has a binomial expression wherein the powers are re- 
placed by the derivatives of the same order. When the function y{x) is constructed up 
to some level of approximation one can, in principle, find cr{x) and substitute for it into 
Eq. (1131. 

We shall consider the implementation of this scheme for four choices of the comparison 
function: 02((t) = l,e2(cr) = a,Q'^{a) = exp(cr) - 1 and e'^{a) = - a^ . At the end 
we shall discuss the relation between the different WKB-type methods and equations of 
the Ermakov-Pinney type. In Appendix O we shall give a pure numerical example of 
application of the method of comparison equations. 

No turning points 

In the case for which the function iJ^ix) does not have zeros (turning points) it is conve- 
nient to choose a comparison function 6^(0") = 1. In this case the function U{a) is simply 
an exponent, while Eq. (14.351) can be rewritten as 

= u.\x)y{x) - ' (4-41) 

dx'^ y \x) 

where the function y{x) is defined by Eq. (14.361) . The above equation is nothing more then 
the well known Pinney or Ermakov-Pinney Eq. [IHlllH], which can be solved perturbatively 
with respect to the parameter e [HI]. We shall give here some general formulae for such 
a solution. It is convenient to rewrite Eq. (I4.4ip as 

a;y - 1 = ey^'y , (4.42) 
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where a "dot" denotes the derivative with respect to x. We shall search for the solution 
of Eq. (lO^ in the form (jiSHD- The zero-order solution of Eq. is 

yo = c^-'/' , (4.43) 
and the general recurrence relation following from Eq. (14.421) for > 1 is 

n— 1 n— 1 



^ n—i n— i ^ ^ 

— ^■■■^i—S[{n-l)-{ki + k2 + k3 + ki)] yk^yk2ykzyk^ 

^0 fei=0 fc4=0 

- 6[n - {ki + k2 + k^ + ki)] yk^yk^yksyki} , (4.44) 



where 6{) is the Kronecker delta symbol. It is now easy to obtain from the general 
expression (I4.44p expressions for particular values of Un 

.1 = 1^ (4.45a) 
= -V^ + h + TT^ (4-45b) 

yl 3^/21/1 , m . 'iym , ^ywi , ^myf /a \ 

= — 2 + + + + ■ 4.45c 

yo yo 4a;2 Auj^yo Au^y^ Au^y^ 

Let us end by noting that the standard WKB approximation corresponds to the case 
discussed here and is, of course, only valid away from turning points. 

One turning point: the Langer solution 

In the case of the presence of a linear zero in the function lj'^{x) one can use the Langer 
solution [37]. In terms of the method of comparison equations it means that one chooses 
the comparison function 0^((t) = a. In this case Eq. (14.351) becomes 

Dividing this equation by (dcr/dx)'^ and differentiating the result with respect to x we get 
an equation which depends only on the derivative & and not on the function a. Such an 
equation can be rewritten in the form 

2uujy^ + Au^yy^ = l + e {3yyy^ + y^^'^y^) , (4.47) 

where the function y{x), as in the preceding section, is defined by Eq. (14.361) . Again we 
shall search for the solution of Eq. (14.471) in the form (14.381) . The equation for yo is 

2d;u;y^, + Aoo^yoy'o = i^'yo) ^ 2/o = 1 • (4-48) 
This equation may be rewritten as 
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and integrated 



3 



\ 2/3 



Finally 



Comparing the coefficients of with n > 1 in Eq. f l4.47p we have 

iuYoVn + {l^uuy'o + 20a;2^ol/o) Vn = 4u;"' {u^YoVn) ' = F^, , (4.52) 
where F„ contains the terms which depend on 1/^,0 < k < n — 1, and can be written as 

n— 1 n— 1 



dxuj (syoVoVo + Vo^Vo) (4.55a) 



Fn = '^■■■'^{-S[n-{ki + k2 + h + h + k5 + ke)]{Auj^yt, + 2ujujyk^)yk^---yk^ 

ki=0 fc6=0 

+ 6[{n-l)-{ki + k2 + k3 + ki + k5 + kQ)] (sjjk.yki + y'ffyk2^ Vk^- ■ -Vke^ ■ (4.53) 
The general expression for ?/„ is 

yn = 4^ / ^^^^n . (4.54) 

Let us also give the explicit expressions for the first terms in the expansion (14.381) 
1 

~ 4^ 

1/2 = J {-^^^^^ylVo - 20t^^?/i?/i?/o - ^Ou'^yoyly^ (4.55b) 

+ ^yimvo + ^yoyiVo + yf\l + ^Vo^yiVo + ^'^ymywi) 

ys = ^^^^ j dxuj {-60uJujy2yiyo - ^Oujuyfy^ - 20Lj^y2yiyQ - 80a;^^ol/2l/i?/o 

-AOu^i/oylyl - AOuj'^yiyly^ - '^Ouj'^ymyo + ^myoyo + Syomyo 
+12yoyoy2yo + ^Vmyo + ^^yiyoyiyl + Uyoymyo + ^^yoyovhl 

+y?yl + ^yf\2yt + ^yf^yiyt + i^y'i^ylyl) • (4.55c) 

Before closing this section we may compare our results with some formulae, obtained by 
Dingle [12]. In the vicinity of the turning point, which we choose as x = 0, the function 
u;^(x) can be represented as 

uj [x) = 7ix + —X + —X + —X H (4.56) 
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while the function a{x) looks like 

a{x) = ao + (Tix + ya;^ + '^x^ + + ■ ■ • (4.57) 

In reference [12] the first five coefficients o"o, ai, o"2, o"3 and (74, as functions of the coefficients 
7i) 72, 73 and 74, were obtained in the lowest approximation by using a system of recurrence 
relations. We can reproduce all these formulae on only using the general formula (14.511) . 
Indeed, on writing 

[ljuj{x)\"' 

using expansions (I4.56P and (I4.57P and comparing the coefficients of different powers of x 
we get the formulae 





= 








1/3 
= ll 






0-2 


1 -2/3 

= ^7i 72 








1 -2/3 
= -^li 73 - 


12 _ 
175^^ 


5/3 2 

72 


(T4 


1 -2/3 

= g7l 74- 


44 _ 
315^^ 


5/3 , 
7273 + 



148 _8/3 3 



2625 



71"^^ 72' , (4.59) 



which coincide with those obtained in . The value of ctq comes directly from Eq. (14.461) 
on ignoring the e-term and using the relations (14.561) and (14.571) . 

We can also find the e-dependent corrections to the function a{x) in the neighbourhood 
of the point x = by using the results (14.591) and the recurrence relations (I4.55al) etc. 
We shall only give here the first correction, proportional to £, to the coefficients do and 
(Ti. The value of ctq comes directly from Eq. (I4.46P using the relations (14.560 . (I4.57P and 

(SSSD 

"° = ' \-Ta 14^ j • ('-'^^ 

In order to obtain ai we introduce into the definition of the function y{x)^ Eq. (I4.36p . the 
expansion (14.571) and 



-1/2 



^2 I ( 3^2CT3 04 ^ , 3 



+ ■■■ (4.61) 



Now, on using the recurrence formula (I4.55ap we find the leading (constant) term in the 
expression for yi 

^ ( ^3 -7/2 , 1 -5/2 CFi(yi^'\ . (Aa'^\ 

1/1 = — I -2^2"^! + l^(^-i(^2(yi I + ■ ■ ■ (4.62) 
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From the formula 

da 1 



(4.63) 



dx (vo + eyiy 

it is easy to find the first correction to the first coefficient ai in the expansion 04.571) 

a, = ^l'^ - 2e^ . (4.64) 

Substituting into Eq. f l4.64p the leading term of yi from Eq. (14.621) . of yo from f l4.6ip and 
using the explicit expressions for the coefficients cTj from Eq. (14.591) we finally obtain 

1/3 I I 3 -11/3 7 _8/3 i _5/3 



(^i = 7i +^(^^^27i ~ 135^37271 +^^47i ' j ■ (4-65) 

The results (14.601) and (I4.65P was found by Dingle [l2], but in his article there were two 
sign errors in the second line of equation (69). 

One turning point: more complicated comparison function 

Let us consider the differential Eq. (I4.32p . when the function u!^{x) has a linear turning 
point, but instead of Langer's comparison function B^(cr) = a consider a more complicated 
comparison function, also having linear turning point, namely 

02((x) = - 1 . (4.66) 
Now, Eq. (14.350 will have the form 

On isolating the function e°", taking its logarithm and differentiating it with respect to 
X we finally get an equation which involves only the derivative &. This equation can be 
written as 

ujY - ^dJooy" - Auo^yy'' + l = e {yy^ - 3yyy^ - y^^^y') . (4.68) 
To find the function y^ one should solve the equation 

uYo - ^dJooy^o - 4c^'?/o2/o + 1 = 0. (4.69) 
On introducing the new variable 

z = LoYo (4.70) 

one can rewrite Eq. fl4.69p as 



iz + l)uj 

^1/2 



(4.71) 
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It is also convenient to introduce the variable 

v = yG . (4.72) 

Eq. fl4.7ip now becomes 

2vH 

which can be integrated, obtaining 



uj (4.73) 



2v — 2arctanw = / cudx . (4.74) 



Finally, for i/q we have the implicit representation 

1 1 /" 

yl arctancui/Q = — / cudx . (4-75) 

Oil 2uj J 

We may now write the equation defining the recurrence relation for yn,n>l 

- AujVoyn+ {^u;Vo - 20^2^01/0 - ^'^'^^vl) Vn = -4 ^""'^°^^"^ = F„ , (4.76) 

where 

n—l n~l 

Fn = ■ ■'^ {S[i^-^)-iki + k2 + h + h)] yk^yk^yk:,yki (4.77) 

A:i=0 fc4=0 

-u;^ S [n-{ki + k2 + k3 + h)] yk^yk^yk^yki] 

n—l n—l 

+ ■ {5 [n-{ki + k2 + k3 + ki + k5 + ke)]{2ujujyk^ + Au"^ yk^)yk^- ■ -yk^ 

A:i=0 fe6=0 

-5 [{n-l)-{ki + k2 + k3 + ki + k^ + k(i)]{?,yk^yk2 + yf^yk^yk^ ■ -^/fce} > 

and 

The solution for ?/„ is 

1 r 

dxuGFn, (4.79) 



Here we give the explicit expressions for the first two functions Fi and F2 

Fi = y^yl-^Moyo-yf^l (4.80a) 

F2 = -GivViyo + ^Ouuyy, + 20uj^ymyt + 40^'2/o2/i2/o + (4.80b) 

+ ^oywl - ^mmyt - 3j/o2/i2/o - ^'^mmyiyl - yf\l - ^yf^yiyt ■ 



62 



Chapter 4. WKB-type approximations 



Two turning points 

Let us now consider the differential Eq. (14.321) . wlien the function u)'^{x) has two turning 
points. This situation, describing the tunnehng through a potential barrier was considered 
in [¥n 1^ . In this case the comparison function can be chosen to be 



e2(a) = a2_a2_ ^4_g^) 

The solutions U{a) of Eq. (I4.33P are the well-known parabolic functions Substi- 
tuting e^{a) of Eq. Km into Eq. (g35|) one can isolate a and after the subsequent 
differentiation one can write down the equation for the function y{x) 



+ a2 - eyy^ = 2djujy^ + Au^yy^ - e {3yyy^ + y^^^y^) . (4.82) 
The lowest-approximation yo can be found from the equation 



2^ujYo + a^ = ^djujy^o + 4^'2/ol/o • (4-83) 



This equation can be integrated and the solution y^ can be written down in the implicit 
form 



o uyo 2 . 
ojyn arcsmh = — udx . (4.84) 



To get the recurrence relations for the functions ?/„,ri > 1 it is convenient to take the 
square of Eq. (IT^ 

AuWy^"^ + 16cj^?/2yio + IQuuj^yy^^ - Au'^y^ - Ac? = 

-e{Ayy^ - Uucuyyy^^ - 2Auj^yfy^ - Adjtoy^^'^ y^^ - Scu^^^^yy^^) 

-e\9fyY + y^'^'y'' + Qy^'^yyy') • (4.85) 



On comparing the terms containing the n*^ power of the small parameter e one obtains 

fyn + gyn = ^ {Gyny=Fn, (4.86) 
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where 

/ = 32uj%y^'^ + IGuu'y'o' (4.87a) 
g = A8u^uVo' + lGOi^Voyo + ^7(^^^'yoyo'-^(^^Vo (4.87b) 

n—1 n—1 

Fn = 4 ^ ■ ■ ■ ^ {uj^6 [n - {ki + k2 + h + ki)] yk^yk2yk3yki 

fcl=0 fc4=0 

-5[{n-l)- {ki + k2 + ks + k^)] yk^ykiykayki} 

n—1 n—1 

- 45^---5^ {5[n-(A;i + --- + M] 

fcl=0 fci2=0 

X [u'^u'^yk^yk^ + iu'^ykii/ki + '^^^^ijk^yki) yk^- ■ ■yki2 
-S[in-l)-{k, + --- + k^2)] 

X (^djuyk^ijk^yk^ + Guj'^ykiyk2yk3 + uju^y^'^^yt^yki + 2uj^yffyk2yk3^ yk^ ■ -l/fcia} 



n-2 n-2 

Y^...J2s[{n-2)-{ki + --- + ki2) 

ki=0 fcl2=0 



X (^yk^yk2yksyki + yf^yt2yk3ykA + ^yt^yk2ykiyk^ yk^- ■ ■yk^2 ■ (4.87c) 



G = exp ( / dx'^ ) . (4.87d) 



/ 

The integration of Eq. (14.861) gives 

1 f , GFn 

yn = -^j dx^ . (4.88) 

The exphcit expressions for the coefficients F„ are rather cumbersome and we shall only 
write down the coefficient Fi 

F, = i2uujmyoyo + '^^^^Mhl + ^^^^yf^l^ + ^uj'^yf^oyo - ^iJoyl ■ (4.89) 



Historical and physical remarks 

We have seen that, for some non-trivial applications of the MCE, instead of the standard 
form of the Ermakov-Pinney equation, one is lead to its generalizations which are rather 
different from the known generalizations of Ermakov systems. Let us note that although 
the Ermakov equation was already reproduced in the MCE context (first of all in the 
important paper by Milne [52]), its generalizations arising in the process of application 
of the MCE [m |12] were not yet studied, at least to the best of our knowledge. We shall 
therefore try to give a short review of the history of the Ermakov-Pinney equation and its 
applications to quantum mechanics and to similar WKB-type problems and discuss the 
physical significance of its non-trivial generalization (for a more detailed account of the 
history of the Eramkov-Pinney equation see e.g. [53]). 
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The relation between the second-order hnear differential equation 



+ uj {x)x{x) = 



(4.90) 



and the non-linear differential equation 

d^p 



a 



p^{x) 



(4.91) 



where a is some constant was noticed by Ermakov [38], who showed that any two solutions 
p and u of the above equations are connected by the formula 



C — 

V? 



a 



(4.92) 



The couple of Eqs. f l4.90p and 04.911) constitute the so called Ermakov system. An impor- 
tant corollary was derived |48j from the formula f l4.92p . Namely, on having a particular 
solution p{x) of Eq. f l4.9ip one can construct the general solution of Eq. f l4.90p which is 
given by 



x(x) = cip{x) exp 




+ C2p{x) exp 




(4.93) 



It is easy to see that the function p{x) plays of the role of an "amplitude" of the function 
x(x), while the integral j[dx/p'^{xy\ represents some kind of "phase". Thus, it is not 
surprising that the Ermakov equation was re-discovered by Milne ^52] in a quantum- 
mechanical context. On introducing the amplitude obeying Eq. fl4.9ip and the phase, 
Milne constructed a formalism for the solution of the Schrodinger equation which was 
equivalent to the MCE. Milne's MCE technique was extensively used for the solution of 
quantum mechanical problems (see e.g. [5^). In the paper by Pinney [lU] the general form 
of the solution of the Ermakov equation was presented, while the most general expression 
for this solution was written down by Lewis [51j . 

A very simple physical example rendering trasparent the derivation of the Ermakov- 
Pinney equation and its solutions was given in the paper by Eliezer and Gray [55J. The 
motion of a two-dimensional oscillator with time-dependent frequency was considered. 
In this case the second-order linear differential Eq. fl4.90p describes the projection of 
the motion of this two-dimensional oscillator on one of its Cartesian coordinates, while 
the Ermakov-Pinney Eq. f l4.9ip describes the evolution of the radial coordinate p. The 
parameter a is nothing more than the square of the conserved angular momentum of 
the two-dimensional oscillator. Thus, the notion of the amplitude and phase acquire in 
this example a simple geometrical and physical meaning. Similarly, the supersymmetric 
generalization of the Ermakov-Pinney equation was considered in the paper by loffe and 
Korsch [561. 
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An additional generalization of the notion of an Ermakov system of equations was 
suggested in the paper by Ray and Reid [SZj- Thy consider the system of two equations 



AH-) 



and 



+ u^{x)p{x) 



1 



/ 



(4.94) 



(4.95) 



where g and / are arbitrary functions of their arguments. The standard Ermakov sys- 
tem fl4.90l) . fl4.9ip corresponds to the choice of functions g{^) = and /(^) = a^. One 
can show that the generalized Ermakov system (14.941) . (14.951) has an invariant (zero total 
time derivative) 



'f,9 



u 



dp du\ 
u- p— 



where 



u 



dx 

f\x)dx 
g{x)dx . 



dx I 



(4.96) 



(4.97a) 
(4.97b) 



The invariant (14.961) establishes the connection between the solutions of Eqs. (I4.94p and 
(I4.95P and sometimes allows one to find the solution of one of these equations provided the 
solution of the other one is known (just as in the case of the standard Ermakov system). 
However, for the case of an arbitrary couple of functions / and g a simple physical or 
geometrical interpretation of the Ermakov system analogous to that given in [55] is not 
known. 

Here we have established the connection between the generalized WKB method (MCE [HI 
1^ ) and the Ermakov - Pinney equation and on the other hand we have obtained a gen- 
eralization of the Ermakov - Pinney equation, different to those studied in the literature. 
One has two second-order differential Eqs. fl4.32p and (14.330 and the solution of one of 
them (14.330 is well known and described. Using the previous convenient analogy [H3] . 
one can say that in this case one has two time-dependent oscillators, evolving with two 
different time-parameters. One can try to find the solution of Eq. (14.321) by representing 
it as a known solution of Eq. (I4.33P multiplied by a correction factor. This correction 
factor plays the role of the prefactor in the standard WKB approach while the known 
solution of Eq. (14.330 represents some kind of generalized phase term. Further, the pref- 
actor is expressed in terms of the derivative between two variables x and a (or in terms 
of the oscillator analogy, between two times). On writing down the equation definining 
this factor, which we have denoted by y{x) (see Eq. (14.361) ) we arrive to Eq. (14.351) which 
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could be called a generalized Ermakov-Pinney equation. For the case when the function 
uj'^{x) does not have turning points the comparison function B^((t) can be chosen constant 
(the second oscillator has a time-independent frequency) and the the equation defining 
the prefactor becomes the standard Ermakov-Pinney Eq. (14.411) . In terms of the two- 
dimensional oscillator analogy [55j this means that we exclude from the equation for the 
radial coordinate the dependence on the angle coordinate a by using its cyclicity, i.e. the 
conservation of the angular moment. 

In cases for which the function uj'^{x) has turning points, as in Sees. 14.4.11 14.4.11 and 
14.4.11 the comparison functions are non-constant and instead of the standard Ermakov- 
Pinney equation we have its non-trivial generalizations. The common feature of these 
generalizations consists in the fact that the corresponding equations for the variable y 
also depend explicitly on the parameter a which does not disappear automatically. To 
get a differential equation for y, one should isolate the parameter a and subsequently 
differentiate ti with respect to x. As a result one gets a differential equation of higher- 
order for the function y{x). Remarkably, for the perturbative solution of these equations 
one can again construct a reasonable iterative procedure. 

Again, it is interesting to look at the generalized Ermakov-Pinney equation as an 
equation describing a two-dimensional physical system in the spirit of the reference [53] . 
In our case one has 



where y plays the role of a radial coordinate, a resembles a phase or an angle (e.g. the 
position of a hand of a clock) and a; is a time parameter. It is important to notice that 
according to the definition of the variable y (14.361) there is a relation between the radial 
and angle coordinates 



which on interpreting a as a phase corresponds to constant (unit) angular momentum. 
In our case however the right-hand side of the radial Eq. (I4.98P contains an explicit 
dependence on the "angle" a. Thus, the couple of Eqs. (14.981) and (I4.99p cannot be 
represented as a system of equations of motion corresponding to some Lagrangian or 
Hamiltonian as in Ref. ^55j. Perhaps, this system could be described in terms of non- 
Hamiltonian dissipative dynamics, but this question requires a further study. 





Chapter 5 



Analysis of Cosmological 
Perturbations 



In this chapter, we shall apply the different WKB approaches reviewed previously to 
the equations that drive the scalar and tensorial cosmological perturbations. In the first 
Section, we recall the main expressions and review one fo the first applications of the WKB 
approximation to cosmolgical perturbations [58]. In the second Section, we present the 
improved version of the WKB analysis [59] and, in the last part, we employ the MCE [60] 

5.1 Standard WKB approximation 

In this Section we begin by recalling that the straightforward application of the standard 
WKB method leads to a poor approximation and that the results can be improved to first 
adiabatic order by means of the transformation (15.91) . as was shown in Ref. [58]. However, 
the extension to higher orders is more subtle, as we shall discuss in some details. 
We have seen that the general equation for cosmological perturbations is 



with j = S for scalar perturbations and j = T for tensorial ones, so Zs = z and zt = a- 
We can now solve Eq. fl5.ip together with the condition that the modes are initially plane 
waves for wavelengths much shorter than the Hubble radius. 




{k,r]) /ij=0. 



(5.1) 



where 



(5.2) 



—i k rj 




(5.3) 
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where d = ds = 1 and d = d^ = y/Sn/ mpi for scalar and tensor perturbations (mpi is the 
Planck mass). So the dimensionless power spectra of scalar and tensor fluctuations are 
given by the usual definitions 



Vc 



2 7r2 



4e 



(5.4) 



The spectral indices and their runnings give the standard parametrization for the power 
spectra. They are usually defined at an arbitrary pivot scale as 



ns - 1 



dlnP, 



d\nk 



dlnVh 



k=k. 



and 



d" Inn 



as 



dink 



d'lnVh 



k=kt 



(dlnfc)^L=,. 
We also define the tensor-to-scalar ratio 

R -- 



(dlnA;)^ 



(5.5) 



(5.6) 



(5.7) 



In the following we shall drop the subscript j, S and T if they are not necessary. On using 
the potentials in Eqs. (13.681) and 03.691) in terms of a and if, one finds, for the quantities 
defined in Eq. 04.51) 

4 

C(e„) ~ sub-horizon scales 
As,T ={ ^ ' (5.8) 



aE 



super-horizon scales. 



The result for super-horizon scales suggests that a change of variable and function are 
necessary in order to obtain better accuracy. 

A convenient change of dependent variable and funtion is the Langer transforma- 
tion 



X = In 



k 



X=(l-ei)^/^e"^'/V, 



(5.9) 



introduced in Ref. [51] is used in order to improve the accuracy. Such a transformation 
bring Eqs. 03. 4p and 03.64p into the form 



dx^ 



+ u;^(x) 



X 







(5.10) 



^Sometimes the spectral indices and their runnings are defined without assuming k = k^. We shall 
always make it clear where we use this assumption. 

^Note that we define x with the opposite sign with respect to that in Ref. f55^. This of course implies 
that our Xi +oo and Xf ^ —oo, whereas Xi —oo and Xf ^ +00 in Ref. 58J. 
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where the (new) frequency u!{x) in general vanishes at the classical "turning point" x = Xq. 
In fact, on employing the transformation (15.91) . one finds (for 6=1) the new equations of 
motion (15.101) . which explicitly read 



dx^ 



dx^ 



+ 



+ 



ei 



2-2ei 



(3-2ei)e2 (l-2ei)e2e3 (l-4ei)e 



ei 



2-2ei 



2(l-ei)2 2(l-ei)3 A{l-e,)^ 

(2 + ei)eie2 



Xs = 
(5.11) 



ei £2 



Ci ^2 



2(l-ei)2 2(l-ei 



4(l-ei) 



XT = . 



From co'g(x) and ti;T(x), respectively given by the expressions in the square brackets of 
Eqs. (15. lip , we now obtain 



rAs,T = e-2-O(e„)~0 



sub-horizon scales 



4 

27 



^s = 7^ ( ei £2 + ei £2 es 



^2 ^3 _|_ ^2 ^3 ^4 



^ + C(e^) super-horizon scales ^^2) 



An 



27 



(eie^ + eie2e3)+0(6^) 



super-horizon scales 



and it is thus obvious that A becomes very small for any inflationary potential satisfying 
the slow-roll conditions [61] in both limits of interest (but, remarkably, not at the turning 
point). As a consequence, the WKB approximation is now also valid on super-horizon 
scales. 

One might naively expect that the transformation (15.91) will improve the WKB ap- 
proach to all orders (in the formal expansion in 6), thus yielding increasingly better results. 
However, since Q diverges at the turning point x = Xq, whether the function /xwkb can 
still be regarded as a solution to Eq. (15. ip for x ~ xq in the adiabatic limit 6—^0 becomes 
a subtle issue. It is precisely for this reason that one usually considers particular solutions 
around the turning point (for example, combinations of Airy functions for a linear turning 
point) which must then be matched with the asymptotic form (14. 2 p on both sides of the 
turning point. The problem with this standard procedure is that there is no general rule 
to determine the actual matching points (which cannot be xq where Q diverges) in such a 
way that the overall error be small. As a consequence, the standard adiabatic expansion 
for 6 <^ 1 may not lead to a reliable approximation around the turning point, which would 
then result in large errors to higher orders and a questionable further extension of all the 
expressions to the interesting case 6 = 1. 

It was in fact shown in Ref. [36] that a simple extension to higher orders of the above 
procedure [22] actually seems to reduce the accuracy with respect to the leading order 
results of Ref. [HE]- We shall see in the next Section that one can overcome this problem by 
introducing mode functions which remain good approximations also around the turning 
point. 
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5.2 Improved WKB approximation 

We consider here the formahsm given in Section 14.31 and apply it to cosmological pertur- 
bations. The basic idea is to use approximate expressions which are valid for all values 
of the coordinate x and then expand around them. It is also easy to show that, for the 
cosmological case, |cr/ci;^| ^ A as given in Eq. fl5.12p for a; ±oo, so that the new 
approximate solutions remain very accurate for large 



5.2.1 Perturbative expansion 

We start by considering the initial conditions and the matching procedure for the pertur- 
bative expansion given in Section 14.31 



Initial and matching conditions 

We first use the initial conditions fl5.3p in order to fix the constant coefficients A±. 
For the function x such conditions become 



lim = d 



2k V 2k 



(5.13) 



where, from the definitions of rj and x, krji ~ —i^i(xi). We then recall the asymptotic 
form of the Bessel functions Ji^ for x — > +oo 



i^iixi 



2 / TT TT \ 

J±mMx)]r^ J^^^^cos[^i{x)T^m--j , 
from which we obtain the asymptotic expression of the mode function in region I, 



(5.14) 



_l_g-*Ci(a;i) 



1 - ei(xi 



2n 



e-"»/2 ^ 



(5.15) 



Ki+f ) - A+ + e~* ) - (A_ + £ A_ J-_ (xo, Xi) +eA+ J++(a;o, Xi)) 
*(3+f)-A+ + e+*(^-^)" {A.+eA.J~_{xo,x,)+eA+J+^{xo,x,))]] . 



The initial conditions (15.131) therefore yield 



A^ = d 



-,2irmT ' 



A_ = -A, 



e J_|__|_(xo, Xij 



k l-Q-iirnn ' " "t 1 + £ J-_ (xq, Xi^ 

Note that A^ does not depend on e, so it holds to all orders. 



(5.16) 



^ Since the initial conditions must be imposed at x in deep region I, we can consider Xi ~ +00. 
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We next impose continuity at the turning point in order to determine the coefficients 
B±. For this we shall need the (asymptotic) expressions of Bessel functions for x 
Xq [Hni [38] 0. In particular, 



„ , , for dC ^ 

r(i±m) ° 



itm 



(5.17) 



rri ±m) 



for X ^ X 



) 



where F is Euler's gamma function. Near the turning point, we just keep the leading term 
in the expansion (14.111) . so that 



. 1 — 2 m I 1 

co'(x)~vC|x — a;o|^" ,,^(x)~2mvC|x — xo|^™ . 
We thus obtain the following forms of the mode functions near the turning point. 



(5.18) 



Xi{x ~ Xq) ~ A+a+{xQ) 



V2m (mVC 



F(l + m) 



(x - Xo) 



V2~rn [mVC 



F(l -m) 



(5.19a) 



Xu{x ^ Xo) ~ fi+6+(xo) 



F (1 + m) 
V2^ (mVC 



(xo - x) 



F(l-m) 

Continuity across Xq then implies that 

i?+6+(xo) = — A+a+(xo) , B_b_{xo) = A_ a_{xo) . 



(5.19b) 



(5.20) 



From the definitions fl4.20p at x = Xq, Eqs. (15.161) and the relations (15.201) we therefore 
obtain 



B_ 



^ ^ ^ e+^'"" + £j++(xo,x,) J 



1 + e J7 ( 



_(xo,Xi) - e J^_(xo,Xi) 



(5.21) 



,+i m 



+ £ J|+(xo,Xi) 1 -ie:^ J__(xo,Xi)/|+(x/,Xo) 



—I e 



\ + e J_,__(xo, X, 
1 + £ t/+_(xo, X. 



1 -«£:/+_ (x/,Xo) 



I) T + 



1 — ie I_^__{xf, Xo) 



/++(x/,Xo) 



"^The Bessel functions are actually regular around x = xq, therefore the following expressions are just 
the leading terms in the Taylor expansion of and I,^. 
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It is apparent from the above derivation that, since the matching between approximate 
solutions in the two regions is always performed at x = Xq, no ambiguity hinders the 
evaluation of the coefficients B± and the accuracy of the approximate solution is therefore 
expected to increase with increasing order. 



Recursive (perturbative) relations 

We can finally obtain expressions for the coefficients a±{x) and b±{x) by solving the 
integral relations (14.201) . We remark that the r.h.s.'s of such equations contain a±{x) and 
bj-{x), which can therefore be determined recursively by expanding them (as well as any 
function of them) in e, that is 

a±{x) = l + J2 ' b±{x) = 1 + J2^' ^± ■ (5-22) 

It will also be useful to define the following integrals 

J7s^'^\xi,X2) = / ^i(^) "-wKy) uis{y) uis{y) dy 



Z^'^\xi, X2) = ^ I au{y)h^^\y)uiis{y)uHs{y)dy , 



where again w, s and s = ±, and g is a non-negative integer. 



(5.23) 



Leading order 



We begin by considering the leading {e^) order 

af{x)=ht\x) = 1 . 
The corresponding solutions are simply given by 



Xi{x) ~ ui{x) = A+ 



UJi(x) 



Xii(x) ~ uii{x) = -A+ 



UJuix) 



{hm [eii(a;)]+e+^^'"-/_„[eii(x)]} 



where we used the zero order forms of Eqs. (15.161) and (15.211) . 

A_ = = -A+ e+^"'" , B+ = -A+ 



(5.24) 

(5.25a) 
(5.25b) 

(5.26) 



We know that these expressions are also good solutions of Eq. (I5.10p near the turning 
point and, deep in region II, we can use the asymptotic form of for x ^ Xf —00 [381 
[39], 

I±m[^u{x)] 



V2vreii(x) 



(5.27) 
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On neglecting the non-leading mode e we then obtain 

(-i _|_ +im.-K\ 



(5.28) 



If we now use Eqs. fl5.4p . fl5.16p and fl5.28p . and the values of we recover the results of 
Ref. [5H] for all quantities of interest, i.e. power spectra, spectral indices and a-runnings 



ns - 1 ~ 3 + 2 



.2?ii,s(fc,r)) 



7reimli\aHj [1 - ei(?7)] u;n,s(^, ^) 
d 6i,s 



V, 



(0) 



d InA; 



as ^ 2 



(d InA:)^ 



and 



16 ( k 



^2 ^u.T{k,ri) 



nml^ \aHj [1 - 61(77)] u;n,T(fc, ??) 



V, 



(0) 



nT ^ 3 + 2 



d In /c 



ttT — 2 



(d InA;)^ 



(5.29a) 
(5.29b) 
(5.29c) 

(5.29d) 
(5.29e) 
(5.29f) 



k=K 



where we have transformed back to the original variables k and r], and all quantities 
are evaluated in the super- horizon limit (i.e. for k -C aH). 

Next-to-leading order 



We now insert the leading order expressions (15.240 into the integral relations (14.200 to 
compute the next-to-leading order expressions 

a^l\x) = J^}_{x,Xi) — e"* 



a«(x) = jflixo, x) -e-'-^^ /+l{xo, x) 



b^_l\x) = -z4"!(x,xo) -2e+''""Ji'!(x,a;o) 



r(o) 



(0) 



(5.30) 



b^}\x) = —iI^L{xf,x) —ie 



(0) 



™/|°|(a;;,a;). 



+{0) 



Jqq =^ J and 1 1 



where we used the simplified notation J, 
also expanded the coefficients A- and B± to zero order in e. 



HO) 
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From Eq. (15.271) . again neglecting the non-leading mode e we obtain the mode 

function at x = xj to first order in e, 



XII (2;/)' 



—zee 

+e 



2 TT UJii{x 



{1 + 



+i m TV 



+ 16 



f) 



/f(x/,a:o) 



.{Xf,Xo) - /|°l(x/,So) 

J\_1{xq, Xi) + J^}^{xq, Xi) I , 

which, using Eqs. (15. 4p . (I5.16p . (15.311) and the values of d, yields all quantities of interest 
to first order in e. For the spectra we obtain 

v,c.vp [i+gf^nxf)] , p.^^r [i+gf^rM] , (5.32) 

where the relative corrections to the leading order expressions (I5.29al) and (15.29dp are 
now given by 

J\_]_{xq, Xi) — J^1{xq, Xi) 

/+! (X, Xo) - 2 if I {x, Xo) + (X, Xq) 



— ee 



jfl{xo,Xi) + /°l{xo,Xi) 



(5.31) 



GREEN/ X 

5^(1) s,T y^) 



+ 



V3 



(5.33) 



S,T 



in which we set e = 1 and n = 1 as required, and write S and T to recall the use of the 
corresponding frequencies. The expressions for the spectral indices and their runnings to 
this order can finally be derived from the definitions (15.51) and (15. 6p . 



Higher orders 

The above procedure can be extended to all orders. The coefficients up to a^f^ and 
b±.^ must be used to compute the integrals J^J^'^^ and I^J^'^^ which determine a±.~^^^ and 
b^^^^\ However, general expressions soon become very involved, and we shall test the 
effectiveness of our method by applying it to several cases of interest. 



5.2.2 Adiabatic expansion 



It is important to observe that the initial conditions (I5.13P and matching conditions at 
the turning point to first order in the adiabatic parameter 6 yield different results with 
respect to the perturbative expansion in e. First of all, the coefficient A^, which was 
left unaffected by the expansion in e [see Eq. (15.161) ]. acquires a correction. The relations 
between A_ , B± and A^ instead remain those given in Eq. (15.260 to all orders, whereas in 
the perturbative expansion they were modified [see Eqs. (I5.16P and (I5.2ip ]. To summarize, 
to first order in 5, one obtains 



A, 



A. 



'''It 



-A+ e 



B 



1-5 



5l(l) 



+ ll{l)[Xi) \ lUJi[Xi 



± 



tA 



± 



(5.34) 
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and the perturbation modes deep in Region II (at x <^ Xq) are finally given by 



Xii(x) ~ Mii(x) {1 + 5 



011(1) (a;) - 711(1) (x) ujii{x) + 



2 Unix) 



-01(1) (xj) +7i(i)(a;i) I - -iuji{xi) 



(5.35) 



We note that, in order to obtain the above result, one must first take the asymptotic 
expansion of the functions u± and then take the derivative, since the reverse order would 
lead to larger errors. The corrected power spectra are then given by 



in which we set 5 = 1 at the end 0. The spectral indices are also given by 



(5.36a) 



ns - 1 = 3 + 2 



d6i,s ^9\{)s 
dink d In /c 



(5.36b) 



riT = 3 + 2 



d ^II,T _^ ^9fv)T 
d In /c d In /c 



and their runnings by 



+ 



d^^fifs 



(dlnfc)^ (dlnA;)^ 



(5.36c) 



^ y(i)T 



Finally, the tensor-to-scalar ratio takes the form 



= 16ei 



.AD 



9(i)T I ^n,s 



(5.36d) 



where all quantities are evaluated in the super-horizon limit {x ^ Xq) and we used the 
results 



fi'(l)S,T(^) 



011(1) (a;) - 7ii(i)(a;) uJu{x) + 



2 Unix] 



01(1) (a^i) 



,(5.37) 



S,T 



where the indices S and T recall the use of the corresponding frequencies. 



^ Since our method has some similarity with that of Ref. [52] (see also Ref. [H]), it is worth pointing 
out that the first order corrections shown here differ from those of Ref. [52] (at least) in that they contain 
contributions from both Regions I and II. 
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Given the complicated expressions for uj, it is usually impossible to carry out the 
double integration that yields analytically, let alone higher order terms. One must 
therefore rely on numerical computations. In the context of power-law inflation we can 
compare results obtained in next-to-leading order from the two expansions with the exact 
result. 



5.3 Method of Comparison Equation 

In order to apply the MCE to cosmological perturbations, we note that the frequency 
uj'^{x) is 

2x 

a;2(x) = ^ - , (5.38) 

with u^lx) given, respectively for scalar and tensor perturbations, by 

Ax) - ' P~^^V I (3-26i)6. {l-2e,)e,e, (l-4e,)ei 

"^^^""^ ~ + 2(l-6i)2 + 2(l-ei)3 +4(l-ei)^ ^ ' 

^^^"^^ " 4 2{l-e,f 2{l-e,r 4(l-ei)^' ^^"^^^^ 

where we omit the dependence on x in the ej for the sake of brevity. From Eqs. (15.381) . 
(I5.39al) and fl5.39bp we can recognise our second-order differential basic equations (15.111) . 
The point x = xo where the frequency vanishes, uj{xo) = (i.e. the classical "turning 
point"), is given by the expression 

xo = \n[u (1-ei)] , (5.40) 

where we have defined u = z/(xo) and ei = ei(xo). We now choose the comparison function 

= - , (5.41) 

and note that then ao = Xq. Solutions to Eq. (14. 6 p can now be expressed, by means of 
Eqs. (I12ZD, (HiSnD and dEUD, as 



X±(x)^ W^J±p(— -) , (5.42) 



uj[x) \ i — ei 



where the J's are Bessel functions [M], and the initial condition (15. 3p can be satisfied 
by taking a linear combination of them. However, in contrast with the improved WKB 
method presented in Section 14.31 MCE solutions need not be matched at the turning 
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point, since the functions (15.421) are valid solutions for the whole range of the variable x. 
Eq. (I4.3ip at the end of inflation, x = Xf, becomes 



-6(^0 - - In 



z/- e(af) 
V + e(af) 



— V 



1 + ln 



1-ei 



- In (2 v) 



(5.43) 



where the super-horizon limit Xf <^ xq (o"f 
One then has 



-oo) has been taken in the second line. 



e'"' 2 V 

— ^ — exp 

1 — ei e 



V 



(5.44) 



Finally, on using Eq. f l5.44p . we obtain the general expressions for the power spectra to 
leading MCE order. 



16 



X=X{ 



eij LOT 



(5.45a) 
(5.45b) 



X=Xl 



where mpi is the Planck mass and the quantities inside the square brackets are evaluated 
in the super-horizon limit, whereas the functions 



[1 - cos (2 TT Z/)] [r(l -I^)] 



2 ' 



(5.45c) 



describe corrections that just depend on quantities evaluated at the turning point and 
represent the main result of the MCE applied to cosmological perturbations |^. The expres- 
sion in Eq. fl5.45cl) is obtained by simply making use of the approximate solutions (15.421) 
and their asymptotic expansion at x cxd to impose the initial conditions (15.31) . In Sec- 
tion l5.2l one finds a similar factor but, in that case, using the Bessel functions of order 1/3 
leads to a large error in the amplitudes. The MCE instead uses Bessel functions of order 
z/, with z/ = 3/2 to leading order in the HFF (i.e. the right index for the de Sitter model), 
which yields a significantly better value for the amplitudes of inflationary spectra. 

The MCE allows one to compute approximate perturbation modes with errors in the 
asymptotic regions (i.e. in the sub- and super-horizon limits) which are comparable with 
those of the standard (or improved) WKB approximation (see Sec. 15. ip . Since these 
methods usually give large errors at the turning point (which produce equally large errors 



^Inside the square brackets we recognize the general resuhs vi^^ and V^'' given by the WKB (leading- 



,(0) 



order) approximation in Section [mi The "correction" g{xo) accounts for the fact that in Eq. (|5.4ip 
is a better approximation than Langer's (MHU]- 
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in the amplitude of the power spectra) it will suffice to estimate the error at the turning 
point in order to show that the MCE is indeed an improvement. To leading order (that is, 
on using the approximate solution fl5.42p ). the MCE gives an error at the turning point of 
the second order in the HFF, which means that we have a small error in the amplitudes 
of the power spectra. Unfortunately, this error remains of second order in the HFF also 
for next-to-leading order in the MCE. We shall see this by applying Dingle's analysis [42j 
for linear frequencies to our case (15.381) . We start by rewriting Eq. (I4.29p as 



- V 



V 


'3 ol 


1 






4 a\ 


~ 2 


^1. 



0, 



(5.46) 



where we dropped the x dependence in uj and a and the order of the derivatives is given by 
their subscripts {p\ = dcr/dx, uf = du'^/dx, etc.). Note that the term in square brackets 
is the error obtained on using the solutions (I5.42p . We then evaluate Eq. (15.460 and its 
subsequent derivatives at the turning point (i.e. at 
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(5.47b) 
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(5.47c) 



f 2(71 e^"^ 1 

(a4(Ti + 3a2a3) 9^ (a) - _ (4 + 18 aa (Xi^ + 7(T3 ai + 12 a^) | 



+ 



87 (T3 a 



(Tn 27 0-40" 



01 '^3'^2 , ^ q'5(T2 , 13 (730-4 1 CTe 

2 (7f 2 (7i 



(7 



(7 



, (5.47d) 



where (a) was defined in Eq. (15.411) and we omit two equations for brevity. In order 
to evaluate the error at the turning point 



3ai 



1 0-3 



4af 2aiJ, 



(5.48) 



we ignore the terms in square brackets and equate to zero the expressions in the curly 
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brackets in Eqs. fl5.47ap - (]5.47dp and so on. This leads to 



a 



0-1 



0-2 



0^3 



ln[(l 



2z/2 

1 



1/3 



5 (2P 



2^1/3 



2\2/3 
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2/3' 



6(2i/3^2y 



3 -21/3 ^2 



+ 



175 (z/2/5c^2)5/3 25 (1^4^2)1/3 175^72 7(21/2^7^2)2/3 



(5.49a) 
(5.49b) 

(5.49c) 
(5.49d) 



and similar expressions for 0-4, cts, and which we again omit for brevity. On inserting 
Eqs. (15381) . fl5:39al) and flKM) in the above expressions, we find the errors to leading 
MCE order 
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TP,S 
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(5.50a) 
(5.50b) 



for scalar and tensor modes respectively. 

On iterating this procedure we can further obtain the errors for the next-to-leading 
MCE order A^p. We first compute next-to- leading solutions to Eqs. f l5.47ap -( l'5.47dll and 
so on by inserting the solutions found to leading order for ai, (T2, . . . , ae into the corrections 
(i.e. the square brackets) and into all terms containing 02 (cr) [12]. This leads to 
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(1) 

TP,S 
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(1) 



TP,T 



31712 

331695 
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ei £2 



€162 



21598 
331695 



£2 £3 



(5.51a) 
(5.51b) 



which show that the next-to-leading MCE solutions lead to an error of second order in 
the HFF, too. We suspect that this remains true for higher MCE orders, since there is no 
a 'priori relation between the MCE and the slow-roll expansions. Let us however point 
out that the above expressions were obtained without performing a slow-roll expansion 
and therefore do not require that the ej be small. 
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Chapter 6 
Applications 



We shall now apply the formalism developed in the previous chapter to some models of 
inflation. We shall start by considering the so called Power-law inflation as a benchmark. 
We shall then proceed by expanding our general expressions to second order in the HFF 
and compare with the results obtained from other approximation methods used in the 
literature. 



6.1 Power- law Inflation 

In this model, the scale factor is given in conformal time by 

a(r/)=£o|^?^+^ (6.1) 

with 13 < —2. The case /3 = — 2 is special, since it corresponds to the de Sitter space-time 
with constant Hubble radius equal to £o = H^^. The frequency is given, both for scalar 
and tensor modes, by 



uj\x) = {l + (3r e'^-^P+-j , (6.2) 

and the horizon flow functions read 

ei = ' e„ = ' n > 1 . (6.3) 

Eq. fIS.ip with the initial conditions 05.31) can be solved analytically and the exact 
power spectra aX Xf ^ —oo are given by (see, e.g. Refs. [HH EHl EHl 166] ) 

(6.4a) 

vrmpj /q 
81 
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where 



m ■■= - 

TT 



r(|/? + i/2|) 

2/3+1 



(6.4b) 



and f{(3 = —2) = 10- The spectral indices and their runnings can also be calculated 
from Eqs. (15. 5p and (15.61) and one finds ns — l = 'n.T = 2/? + 4 and as = ot = 0. Exact 
scale-invariance is obtained for /3 = —2. For this particular case, we plot the frequency in 
Fig. 16.11 the function ^{x) in Fig. 16.21 and the perturbation a{x) in Fig. 16.31 




Figure 6.1: The frequency u"^ for (3 = —2. The turning point is at x=k = ln(3/2) ~ 0.4 




^The case /3 — —2 is singular (ei = and the expression of the scalar power spectrum blows up) and 
should be considered separately. However, there are no scalar perturbations in De Sitter space-time to 
first order. 
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Figure 6.3: The "perturbation" a for /3 = —2. 
6.1.1 Improved WKB leading order results 

The leading order expressions in the perturbative expansion of Section 14.3.11 as well as 
the adiabatic expansion of Section I5.2.2[ lead to the same result as in Ref . [58j , 



j<°>(/J) *: 



2/3+4 



(6.5a) 



2 62/3+1 



|2/3+2 • 



where the function g^'^\P) is given by 

\2/3 + l\ 

This yields a relative error for the amplitude of the power spectrum 



.(0) .. 



100 



/(/?) 



% 



(6.5b) 



(6.6) 



which decreases for increasing \j3\, as can be seen from Fig. 16.41 but is rather large (about 
10%) for the de Sitter space-time. The spectral indices and their runnings are instead 
predicted exactly as ris — 1 = nx = 2 /3 + 4 and as = = 0. 



6.1.2 Improved WKB next-to-leading order results 

We shall now compare the corrections to the amplitude of the power spectra coming from 
the two different expansions we described in Chapter O 

Perturbative expansion 

From the next-to-leading expressions in Section 14.3.11 we obtain that the relative cor- 
rection to the power spectra is given by the function (^^^^^^ := = fi'(i)^T™ 
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Ap 




-10 -8 -6 -4 p -2 



Figure 6.4: The total error on the power spectrum to leading order and to next-to-leading 
order for the perturbative correction (superposed upper lines) and to next-to-leading order 
for the adiabatic correction (lower line). 

Eq. fl5.33p evaluated ai x = x j x^:. This quantity can be determined numerically and 
one finds that it approaches an asymptotical finite value for xj — oo. Some examples 
are given in Table EH] for Xf = —13, from which it is clear that the improvement over the 
leading order is small. 
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+7.1 ■ 10^^ 
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3.6% 
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2.5% 


^(1) 


10.4% 


6.3% 


4.6% 


3.6% 


2.9% 


2.5% 



Table 6.1: Next-to-leading order improvement for the power spectra (fi'(^)^^^'~^) and total 
final error (Ap'') to first order in e. The latter is essentially the same as to leading order 



Once the correction has been included, the total relative error on the power spectra 
(denoted by Ap'*) therefore remains of the same order as the leading order error Ap'*. 

Adiabatic expansion 

We analogously evaluate the coefficients and 7(i) numerically. For the case P = —2, 
the result is plotted in Fig. 16. 5[ from which it appears that the functions 0(i)(x) and 
7(1) (x) converge to finite asymptotic values for x — > ±oo. Since the same behavior is 
also found for all (3 < —2, we exhibit the values of such coefficients (with Xj = 20 and 
Xf = —13) in Table [6^ together with the corresponding next-to-leading order relative 
improvement [g^{xf) as defined in Eq. (15.371) ] and the total error on the power spectra 

On comparing with the last row in Table 16.11 we can therefore conclude that the adia- 
batic corrections are significantly better than those obtained from the perturbative expan- 
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Figure 6.5: The functions (p(i){x) (solid line) and 7(i)(x) (dotted line) for /? = —2. 



sion and, in fact, yield the amplitudes with extremely high accuracy (see also Fig. 16.41) . In 
fact, this method seems (at least) as accurate as the next-to-leading order in the uniform 
approximation employed in Ref. 





-2 


-3 


-4 


-5 


-6 


-7 


01(1) 


+3.2-10-3 


+ 1.1 ■ 10-3 


+6.0-10-^ 


+3.6 ■ 10-^ 


+2.4- 10-^ 


+1.7-10-^ 


7i(i) 


-2.4 ■ 10-19 


-5.9 ■ 10-20 


-2.6 ■ 10-20 


-1.5-10-20 


-9.4 ■ 10-21 


-6.6 ■ 10-21 


011(1) 


+4.5 ■ 10-3 


+1.6 ■ 10-3 


+8.2 ■ 10-^ 


+5.0 ■ 10"^ 


+3.3 ■ 10-^ 


+2.4 ■ 10-^ 


7ii(i) 


-3.5 ■ 10-2 


-1.2 ■ 10-2 


-6.4 ■ 10-3 


-3.8 ■ 10-3 


-2.6 ■ 10-3 


-1.8 ■ 10-3 


AD 
9(1) 


+1.1 • 10-1 


+6.3-10-2 


+4.5 • 10-2 


+3.5 ■ 10-2 


+2.9 ■ 10-2 


+2.4 ■ 10-2 




0.83% 


0.50% 


0.35% 


0.27% 


0.22% 


0.18% 



Table 6.2: First order coefficients 0(i) and 7(i), correction 
(Ap'*) for the power spectrum to first order in 6. 
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6.1.3 MCE results 



The MCE (developed in Scection 15.31) yields the exact power spectra, spectral indices and 
runnings. 
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where ipi = nip^ is the Planck length and 
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(6.7c) 
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with r the Gamma function. The spectral indices are ns — 1 = t^t = 2/3 + 4 and their 
runnings ctg = ax = 0. Finally, the tensor-to-scalar ratio becomes 

fi=16|±|, (6.8) 

which is constant as well. 



6.2 Slow- Roll Inflation 



Let us briefly recall the known results obtained from the slow-roll approximation [671 EHl 
. To first order in the slow- roll parameters fl3.65l) . we have the PS 



re ei mpj 



1-2(C+I)ei-Ce2- (261 + 62) In ( ^ 



(6.9a) 



77 m 
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1 -2(C+ l)ei -2 61 In ( — 



where C = 7g+ln2 — 2^ —0.7296, 7^, being the Euler-Mascheroni constant, and all 
quantities are evaluated at the Hubble crossing 77=,,, that is the moment of time when = 
{aH){N^) = {aH)^:. (The number is usually called "pivot scale".) From Eqs. fl6.9al) . 
we can also obtain the spectral indices and a-runnings. 



ns-1 



-2ei - 62 , riT = -2 61 



(6.9b) 



as = Qt = 



(6.9c) 



on using respectively Eqs. (15.51) and (15.61) . From Eq. (15. 7p . the tensor-to-scalar ratio 
becomes 



i?= I661 



1 + C e2 + 62 In 



k 



(6.9d) 



6.2.1 Improved WKB leading and second slow- roll order 

We now consider the possibility of obtaining consistent second order results in the param- 
eters 6j's from the leading WKB approximation. 

We first set g^^ix) = in Eqs. flCTil) . fl536bl) . flHIM and fl536dl) . We then find it 
convenient to re-express all relevant quantities in terms of the conformal time t]. For this 
purpose, we employ the relation 



kr] 



aH{l - ti 



[1 + 6162 + 0(6^)] 



;i-ei) 



[1+6162 + 0(6^] 



(6.10a) 
(6.10b) 
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We can now expand the frequencies to second order in the HFF, 

/g 3 \ \ 



'9 1 
^^t(^) = fc'r/' (l-2eie2)- ( - + 3ei+4e2_-e,e2 



and the PS become 

^ i/? 2 / 8 1 19 2 19 



(6.11) 



(6.12) 



_ 16i/f2 , ,3 2 / 8 19 . 26 , 



From now on, the subscript / will denote that the given quantity is evaluated in the 
super- horizon limit. 

We must now compute the arguments of the exponentials in the above expressions, 
which are in general of the form 



no 



dry 



^n{m,Vuk)= ^A^ir^)-k^r,^-L , (6.13) 

Jrjf V 

where the function ^^(77) contains the HFF, but does not depend on k, and we have used 

dx = ^ [l + e,{ri) e^iv) + 0(ef)] . (6.14) 

Let us also note that, at the turning point rj = rjo, one has 

A{7]o) = -k7]o. (6.15) 

It is now clear that, in order to obtain consistent results to second order in the slow-roll 
parameters, we must consider the time-dependence of the e^'s, and the function A'^{ri) may 
not be approximated by a constant. This does not allow us to perform the integral unless 
the scale factor a = a{ri) is given explicitly or, as we shall see, some further approximation 
is employed. 

In Refs. [TOl [71] , we proposed a procedure which will now be described in detail. Let 
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us start with the general exact relation 
m 
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A{r]) - ^A^{r])-k^r]A [A\r])]' 
A{r]) + y^{r]) - k^ 7]^ ) 4:A{r]) 



dr] , (6.16) 



which holds for every function A'^{t]). We can derive this relation by initially considering 
the integral in Eq. (16.131) with A constant (i.e. independent of 1], which is just power-law 
inflation). In this case, the integration can be performed and yields the first and second 
term in the r.h.s. of Eq. fl6.16p with the function A independent of t], whereas the third 
term vanishes identically. If we then reinstate the r^-dependence of A in the two non- 
vanishing terms and differentiate them with respect to r/, we obtain the original integrand 
in Eq. (I6.13P plus another term, originating from A'^{r]). We then obtain the result (16.161) 
by (formally) integrating the latter term and subtracting it from the previous two. 

One can also repeat the above procedure for the new integral in the r.h.s. of Eq. (16.161) . 
To do this, we define 



[A\r])]'-: 

and, after several simplifications, we can write the general relation as 
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and Li2(-2) is the dilogarithm function (see, e.g. Ref. 

z'' _ r In f 1 - z') 

k=l 



Jo ^ 



On putting together all the pieces, we finally obtain 



+ ^ln 
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dr] . 



For the cases of interest to us [i.e. for the frequencies (16. lip ], the functions A for scalar 
and tensor perturbations to second order in the HFF are given by 
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(6.22a) 
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+ 3ei + 4e^ + 4eie2 



so that 



-Bs = 3 ei 62 + - 62 €3 , -Bt = 3 ei €2 



(6.22b) 



For the results to second order in the HFF, we can neglect the last two integrals in 
Eq. fl6.2ip . since B'i^r]) = O (ef) and B'^{ri) = O (ef). Thus our approximation consists in 
neglecting such terms and not in assuming the HFF be constant (or any specific functional 
form for them). We can finally write ^ii,s('7o5 Vfi ^) C,ii^t{vo, Vfi ^) 



3 3 1n3 , ln3 

4ii,s^- 2 2 ^ ^^'^ ~ ^^'^ 
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3 3 1n3 , / 1\ 2 /4 1n3 3 tt^A 
6i.T— - + ^ + ln3e,, + ( ln3 + - 1 e,, + i— - + - e,,e,, 



+ - ei,f - + (^In 3 - ^ ) ei,f eaj 
ei_f e2,f In^ (-/cr^f) , 



In {—krif) 



(6.23b) 



where the arguments in have been omitted. On inserting the above quantities in 
Eqs. f l6.12p . we can calculate the PS. 

In order to compare the PS from Eqs. 06.121) . fl6.23ap and fl6.23bp . with the slow- roll 
expressions in Eqs. f l6.9al) . we need a relation between H{N{) and if(A^*), ei(A^f) and 
ei(A^=K), and so on, to second order in the HFF (hereafter, quantities without a subscript 
will be evaluated at the Hubble crossing r^* corresponding to the pivot scale K). We 
expand the parameters e^^f to first order in AA^ = Nj — N^, in the numerators and to 
second order in the denominator of the scalar spectrum. 



-) 1 

^ ~ 1 + e,+i AiV + - {el, + e,+i e.+s) AiV^ 



(6.24) 



and, analogously, we find 



1 - 2 ei AAT - (ei €2 - 2 e^) AA^^ . 



(6.25) 



We can eliminate rji in the logarithms by expressing it in terms of l/{aH)i 
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(6.26) 



where in the first equality we used Eq. (I6.10al) . in the second one the definition of the 
pivot scale, and, in the last relation, Eq. (16.241) . If we now use the above expressions we 
obtain the following PS, where the superscript stands for second slow-roll order and 
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the subscript WKB for leading adiabatic order, 



V, 



(2) 

CjWKB 



A. 



1 - 2 (DwKB + 1) ei - Z^wKB e2 + ( ^ I^wKB + ^ ) ^2 



20 



+ -DwKB + 2 -DwKB - - J ei + ^-DwKB ~ -DwKB ^ ~ ) ^2 

+ ^-2^wKB + ^-Y^J e^es 

+ [-2ei-e2 + 2(2DwKB + l) + (2 Z^wkb - 1) £162 
+-DwKB £2 - -CwKB £2 £3] In {^Y^ 

+\ {Ul + 2e,e, + el-e,e,) In^ } 



(6.27a) 



V, 



(2) 

/l.WKB 



A. 



TT m 



PI 



l-2(DwKB + l) ei+{2D 
19 \ 



2 

WKB 



2 -DwKR 



2 A 



+ 



12 



+ [-2ei + 2(2DwKB + l) e?-2 (Z^, 
In the above, we have defined -Dwkb = — ln3 



l).,..]ln(±) 



I — In 3 ~ —0.7653, which approximates the 
coefficient C in Eqs. f l6.9ap with an error of about 5%, and the factor A^kb = 18/e^ ~ 
0.896 which gives an error of about 10% on the amphtudes. An important resuh is 
that the PS to second order in the slow-roll parameters do not depend on AA^, once the 
laborious dependence on the Hubble crossing has been worked out: this fact is in complete 
agreement with the constancy of the growing modes of ( and h on large scales. The 
spectral indices, from Eqs. (15.36b[) . and their runnings, from Eqs. fl5.36cp . are analogously 
given by 
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From Eq. (15.36dp the tensor-to-scalar ratio becomes 
-RwKB = 16 ei 



1 + -DwKB ^2 + ( -DwKB + Cl ^2 + -D^KB ~ ^ 

1 „o TT^ 1^ 



24 1^ 



+ (e2 + ei 62 + -DwKB 4 + -DwKB £2 63) In 

4 ('^^--)>"'(^)^ 
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We note that Eqs. fl6.27al) . f l6.27bl) and f l6.27cp agree with Ref. [58], to first order in the 
HFF. 

We can also compare the above results with others as we show in Appendix [Dl For 
example, one could approximate the HFF by their Taylor expansions around rjo as in 
Ref. [731 ESI Clj, and give the results in terms of the e^.'s calculated at that point (see 
Table [D^ for results only on spectral indices). 

6.2.2 Improved WKB next-to-leading and first slow-roll order 

To first order in the HFF, the function A can be taken as constant. In this case, let us 
then introduce the more convenient notation 



$,±{xo,x;k) = ± / \/T^^^iy)dy 

J X 



= ±A [atan± (e±) - 0±] , (6.28) 
where we have used dx = dr/r, and 
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^U^i) = T + 3ei + ^e2 



(6.29a) 



krio = -A, (6.29b) 
with ei and €2 constant as well. It is also useful to define 



e±^^/±(l-^) (6.30a) 
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atan+(x) = arctanh(a;) 

(6.30b) 

atan_(a:;) = arctan(a;) , 

where the plus (minus) sign corresponds to region II (I). With this notation, the expres- 
sions for 7(i)(?7) and 0(i)(r/) in Eqs. fl4.25p can be explicitly written as 



1 r 5^30 



7±(i)(^) = -T^ <^ ^7^S^± 



181 



where we have used 
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The limits of interest are then given by 

lim 7_(i)(r7) = (6.32a) 

23 

^-<.,('rt = (6.32b) 

lim 7+(i)(^) = - TTT-Ti (6.32c) 



lim 0+(i)(r7) = ^^^37^ , (6.32d) 



lim 6_(r7) = cx) (6.32e) 

?7— > — oo 



lim 6+(r/) = 1 , (6.32f) 



lim (uju + ^] =A. (6.32g) 

We can now use the values in Eqs. f l6.29al) to calculate the correction in the super- 
horizon limit [r] —>■ 0^) for scalar and tensor perturbations, obtaining 

AD _ 37 19 / 1 
" 324 " 243 



(6.33) 
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Finally, from Eqs. fl5.36al) and fl6.33p . we can write the expressions for the scalar and 
tensor PS, with the superscript for first slow-roll order and the subscript WKB* for 
next-to-leading adiabatic order, as 



V, 



(1) 

(,WKB* 



A. 
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(6.34a) 



where -D^kb. = ^ ~ ^^^^ ~ —0.7302, which approximates the coefficient C in Eqs. f l6.9al) 
with an error of about 0.08%, and the factor A^kb* = 361/18 ~ 0.999 which gives an 
error of about 0.1% on the amplitudes. We also obtain the standard slow-roll spectral 
indices and a-runnings, 
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on respectively using Eqs. (15.36bl) and fl5.36cp . From Eq. ( ]5.36dl) the tensor-to-scalar ratio 
becomes 
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1 + ^wKB. £2 + £2 In ( 7^ 



(6.34d) 



6.2.3 Improved WKB next-to-leading and second slow-roll or- 
der 



In order to give the results to next-to-leading WKB order and second slow-roll order, we 
should evaluate the corrections for scalar and tensor perturbations to second order 
in the e^'s. This, in turn, would require the computation of the functions 0±(i) and 7±(i), 
in Eqs. (14.251) . whose analytical calculation is extremely difficult. We shall therefore try 
to give the expressions for the corrections by following a heuristic (and faster) method. 

Let us start from the first order slow-roll expressions given in Eqs. fl6.33p and add all 
possible second order terms in the HFF. Some of such terms can then be (partially) fixed 
by making reasonable requirements, which will be explained below, so that the corrections 
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First of all, we have factorized the number 361/324 so as to recover the standard slow- 
roll amplitudes fl6.9ap with a very good accuracy. As already mentioned, the first order 
terms are the same as those in the results f l6.33p of Section 16.2.21 evaluated however 
in the super-horizon limit. The form of the coefficients multiplying the second order 
monomials in the HFF have been partially determined by using the expressions given 
in Section ESSl [Eq. ^M^-^^^] with Eqs. ^:TM replacing Eqs. ^^IM) . By this 



procedure, one expects to obtain results correct up to derivatives of A , which can just 
contain mixed terms to second order. The numerical coefficients in front of e\ j and e\ f 
are thus uniquely determined, whereas the coefficients multiplying the mixed terms ei^f e2,f 
and e2,f ea^f remain partially ambiguous. The last number in each square bracket arises 
directly from this procedure, whereas the other terms are chosen so as to match the 
"standard" dependence on Z^wkb. in the PS, with 6s,t and rfg left undetermined. Finally, 
we have also required that the PS do not depend on AA^, consistently with the constancy 
in time of the growing modes for h and C,. This requirement fixes the coefficients in front 
of the logarithms in Eqs. fl6.35p uniquely 0, and leads to spectral indices which do not 
depend on 6s,t and do,. 

Proceeding then as in Section 16.2.11 from the above corrections we obtain the expres- 



^For the same reason (i.e. to avoid the appearance of AiV), we cannot consider terms such as 
O (ef) In^ i-krii) in Eqs. 
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sions for the scalar and tensor PS, 
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We also obtain the spectral indices fl5.36bp and their 



runnings fl5.36cp . 
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and the tensor-to-scalar ratio becomes 
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Let us end this Section with a few remarks. The undetermined coefficients 6s,t and 
ds still appear in the PS (and their ratio R), but not in the spectral indices and runnings, 
which are therefore uniquely specified by our (heuristic) procedure. In general, one may 
expect that the complete adiabatic corrections to second order contain HFF calculated 
both in the super-horizon limit and at the turning points (zeros) of the frequencies. How- 
ever, this is not relevant in the present context, since the possible numerical coefficients 
that would multiply such terms can be eventually re-absorbed in the definitions of 6s^t 
and (ig. This could only be an issue (in particular for the spectral indices) if we considered 
third order terms (since it would then become important where the second order terms 
are evaluated). 

6.2.4 MCE and second slow-roll order 

We now consider the results fl5.45ap - fl5.45cl) given by the MCE to leading order (denoted 
by the subscript MCE) and evaluate them to second order in the HFF (labelled by the 
superscript (2)) for a general inflationary scale factor. A crucial point in our method 
is the computation of the function ^ defined in Eq. f l4.3ip which can be found in detail 
in Section HI of Ref. [71]. For the sake of brevity, we shall not reproduce that analysis 
here but it is important to stress that, in contrast with the GFM and other slow- roll 
approximations, it does not require a priori any expansion in the HFF since Eq. (34) of 
Ref. [T^ is exact and higher order terms are discarded a fortiori. From that expression, 
upon neglecting terms of order higher than two in the HFF, we obtain the power spectra 



a. = ^{l-2(C+l).,-C..+ (2C^ + 2c4- 




+ [-2ei-e2 + 2(2C + l)e? + (4C-2DMCE-l)eie2+Ce2-D^cEe2e3] In 





(6.37a) 
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where Dmce = | - ln3 ^ -0.7652 and C = In 2 + 7e - 2 -0.7296, with 7e the Euler- 
Mascheroni constant. A clarification concerning the g{xo)^s is in order. Since the turning 
point does not coincide with the horizon crossing where the spectra are evaluated [7l] , we 
have used the relation 

eiixo) c^ei-ei e^+i In , (6.38) 

in order to express the g{xo)^s as functions of the crossing time (HFF with no explicit 
argument are evaluated at the horizon crossing). The spectral indices are given by 

"s!Le - 1 = -2 ei - 62 - 2 Ci - (2 L'mce + 3) ci £2 - /^mce £2 £3 (6.39a) 

4'lcE = -2 ei - 2 - 2 (Dmce + 1) d £2 , (6.39b) 
and their runnings by 

(2) 

"s,MCE = -2 ei £2 - 62 63 (6.39c) 

«T,LcE = -2 ei £2 • (6.39d) 
The tensor-to-scalar ratio becomes 

ff - 1-^C^=-h(c^-tH-5).....(1c^-^ + i).? (6.40) 

C Z)„c, - i Dl„ + C hi(2) - £>„c, lri(2) - ^ + i) . 

The polynomial structure in the HFF of the results agrees with that given by the Green 
Function Method (GFM) [7S1 P] and the WKB approximation [Ml EHl EHl EI] (the same 
polynomial structure is also found for the spectral indices by means of the uniform ap- 
proximation pi]). Let us also note other aspects of our results: first of all, the factors 
g{xo) modify the standard WKB leading order amplitudes [Ml EH] so as to reproduce the 
standard first order slow-roll results; secondly, we found that C and -Dmce are "mixed" in 
the numerical factors in front of second order terms (we recall that -Dmce differs from C 
by about 5 %); further, -Dmce = -Dwkb of Refs. [Ml IMl EHl EI]- The runnings as and 
are predicted to be C(e^) [76], and in agreement with those obtained by the GFM [751168]. 



6.2.5 Second slow-roll order MCE and WKB versus GFM 



We shall now compare the slow-roll results obtained from the MCE and other methods 
of approximation previously employed [71] with those obtained using the GFM in the 
slow-roll expansion. We use the GFM just as a reference for the purpose of comparing 
the other methods among each other and of illustrating deviations from pure slow-roll 
results. For a given inflationary "observable" Y evaluated with the method X, we denote 
the percentage difference with respect to its value given by the GFM as 



100 



Y - Y 

X GFM 



(6.41) 



6.2. Slow-Roll InHation 



99 



where X = WKB stands for the first WKB and second slow-roll orders [70j, X = WKB* 
for the second WKB and second slow-roll orders [7T] and, of course X = MCE for the 
result obtained in this paper. Note that for the case X = WKB* we shall set the three 
undetermined parameters 6s = &t = c^s = 2 in order to minimize the difference with 
respect to the results of the GFM (see Appendix |E] and Ref. [7lJ). In Fig. 16.61 we show, 




A f? 



Figure 6.6: Percentage differences (16.411) between scalar-to-tensor ratios given by the GFM 
and those obtained from the WKB, WKB* and MCE. 

with dots of variable size, the percentage differences (16.411) for the scalar-to-tensor ratios 
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R at the pivot scale A; = /c* for < ei < 0.5, |e2| and {e^l < 0.5. From the plots it appears 
that the level of accuracy of the MCE is comparable to that of the WKB* and both are 
(almost everywhere) more accurate than the WKB. However, the MCE achieves such a 
precision at leading order and is thus significantly more effective than the WKB*. In 
Fig. 16.71 we show the difference in ns and the relative difference in P(^{k^) between the 
MCE and GFM. 




I X 100 APcmce(A;*) 



Figure 6.7: Left box: absolute difference between scalar spectral indeces ng evaluated 
with the MCE and GFM (rescaled by a factor of 100 for convenience). Note that relevant 
differences of order 0.05 (shown as 5) occur at the boundaries of the intervals considered for 
the ej's. Right box: percentage difference 06.411) for the amplitudes of scalar perturbations 
at the pivot scale = fc* between the MCE and GFM. See Fig. l6.6l for the meaning of 
dot size. 

In Fig. 16.81 we finally plot the relative difference in Pc_{k^) for 62 = — 2ei (the scale 
invariant case to first order in slow- roll parameters). 

6.2.6 GFM and slow-roll approximation 

Let us further compare our results with those obtained by the GFM. As we stated before, 
the most general result of the present work is given by the expressions for the power 
spectra f l5.45al) and ( ]5.45b[) with the corrections f l5.45cp . In fact, the spectral indices 
and their runnings are the same as found with the standard WKB leading order [581 E] • 
On the other hand, the main results of the GFM are the expressions for power spectra, 
spectral indices and runnings ". . . in the slow-roll expansion." that is, for small HFF 
(see, e.g., Eqs. (41) and (43) in Ref. [75]). Here, and in some previous work [58], [59] . 
we instead obtain results which hold for a general inflationary context, independently of 
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Figure 6.8: Percentage difference (16.411) between the MCE and the GFM for P^{K) re- 
stricted on the hypersurface €2 = — 2ei, which corresponds to a scale-invariant spec- 
trum to first order in the slow-roll expansion. The graph is given for < ei < 0.5 and 
-0.5 < £3 < 0.5. 



the "slow-roll conditions" |f| or "slow-roll approximation" B The slow-roll case is just a 
possible application of our general expressions which can be evaluated for any model by 
simply specifying the scale factor. Our general, and non-local, expressions in fact take 
into account all the "history" of the HFF during inflation. For example, we note that the 
MCE at leading order reproduces the exact result for power-law inflation (see Sec. 16.11) . 
whereas the GFM reproduces those to the next-to- leading order [75] . 

Chaotic inflation [77] is another example for which a Taylor approximation of the HFF 
(such as the one required by the GFM [75] or our Green's function perturbative expansion 
in Appendix [Fj may lead to inaccurate results. We consider a quadratic potential 



V{<j>) = -m'<j>' 



(6.42) 



where m is the mass of the infiaton 0. In a spatially fiat Robertson- Walker background, the 
potential energy dominates during the slow-rollover and the Friedman equation becomes 



An m 



2 a2 



(6.43) 



The Hubble parameter evolves as 



H 



An V 



m 



PI 



(6.44) 



On using the equation of motion for the scalar field 



(6.45) 



■^Our approximation method does not require small HFF. 
We expand instead according to the scheme given in Ref . [TT] . 
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6- 

H 

4^ 



5 10 15 20 ^ 25 

Figure 6.9: Numerical evolution of H (solid line) and its analytic approximation (dashed 
line). The initial condition corresponds to Hi ^ 8.2m (time in units of 1/m) when the 
inflaton starts in the slow-roll regime with a value 4mpi. 

and neglecting the second derivative with respect to the cosmic time, we have 3if0 ~ 
— 0. Eq. (16.441) then yields 

2 

TTl 

H^-—=Ho, (6.46) 

which leads to a linearly decreasing Hubble parameter and, correspondingly, to an evolu- 
tion for the scale factor which is not exponentially linear in time, i.e. 



H{t) Ho + Hot , (6.47a) 
a{t) ~ a(0) exp ( Ho t + Ho^-] = a(0) exp [N{t)] . (6.47b) 



In Fig. l6.9l the analytic approximation fl6.47al) is shown to be very good on comparing with 
the exact (numerical) evolution. We are interested in the HFF in the slow-roll regime. 
From the definitions f l3.65p we then have 

^iW = ^^-^ (6-48a) 

Ho + Hot 



e-W = ,2 = 2 ei(t) , n > 2 , (6.48b) 

Ho + Hot 



which we plot in Fig. I6.1U[ 

Let us take the end of inflation at tend — 24/m (so that the analytic approxima- 
tion fl6.47ap . (16.47bp is very good till the end), corresponding to A^'end = 125 e-folds. We 
then consider the modes that leave the horizon 60 e-folds before the end of inflation and 
take that time {t = t^ ~ 8/m or iV* = 125 — 60 = 65) as the starting point for a Taylor 
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Figure 6.10: HFF in the slow-roll regime: the dashed line is ei and the solid line is all the 
e^'s with n>2 (time in units of 1/m). 



expansion to approximate the HFF 0, 



ei (AN) ~ + Ci+i AN+-ei (e^+i + e^+i e^+a) AA^^ 



(6.49) 



where the e^'s in the r.h.s. are all evaluated at the time t^, and AA^ = N — A^^,. The 
percentage error with respect to the analytic expressions (I6.48ap and (16.48bp , 



& = 100 



e,(t) - (AN) 



(6.50) 



with AA^ = N(t) — N^: is then plotted in Fig. 16.111 for i = 1. It obviously becomes large 
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Figure 6.11: Percentage error for ei to first and second order in AA^ (dashed and solid 
line respectively; time in units of 1/m). 

rather quickly away from t = t^, and we can immediately conclude that, had we used 
a Taylor expansion to approximate the HFF over the whole range of chaotic infiation 



^One can of course conceive diverse expansions, e.g. using other variables instead of N ^ however the 
conclusion would not change as long as equivalent orders are compared. 
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(—65 < AA^ < 60), we would have obtained large errors from the regions both near the 
beginning and the end of inflation. 

In general, we expect that the Taylor expansion of the HFF will lead to a poor de- 
termination of the numerical coefficients (depending on C in Eqs. f l6.37al) and ( ]6.37b[) ) 
in the second slow-roll order terms for those models in which the HFF vary significantly 
in time. In fact, we can calculate the integral of Eq. (21) in Ref. [75] with ei instead of 
the complete ^'(Inx) and yo{u) instead of y{u). The percentage difference between this 
integral calculated using the Taylor expansion (16.491) with respect to the same integral 
calculated with ei in Eqs. (I6.48ap and ( ]6.48bl) is 92 % and 88 % to first and second orders 
in AA^, respectively. The MCE method does not require such an expansion and does 
therefore not suffer from this restriction. 



. 92 
. 9 

Log(PQ) 0-88 




Figure 6.12: Spectrum of the Mukhanov variable Q {Pq = |Qfcp/(2 7r^)) evaluated at 
the end of inflation for the chaotic model of Eq. (I6.42p . The dots represent the numerical 
values and the solid line the analytic flt based on Eqs. (I6.39al) and (I6.39cl) obtained by 
the second slow-roll order MCE approximation. The flrst order and GFM analytic results 
are not shown since they are almost indistinguishable from the MCE result plotted. K 
crosses the Hubble radius at 0=,, ~ 3mpi (the lines are normalized at 10^'^^ A;=k). 

We have also compared our analytic MCE results in Eqs. (I6.39ap and fl6.39cp with 
the numerical evaluation of the spectrum for scalar perturbations. In particular, we 
have evolved in cosmic time the modulus of the Mukhanov variable Q (recalling that 
Pc = |(5a:P/(2 TT^)), which satisfles the associated non-linear Pinney equation [78]. The 
initial conditions for the numerical evolution are flxed for wave-lengths well within the 
Hubble radius and correspond to the adiabatic vacuum, i.e. \QkiU)\ = ^/{o-iU) Vk) and 
\Qk\(ti) = —H(ti) \Qk(ti)\ [79]. The agreement of the analytic MCE approximation with 
the numerical results is good, as shown in Fig. 16.121 

6.2.7 Beyond slow- roll 

The slow-roll approximation is quite accurate for a wide class of potentials. However, 
violations of the slow-roll approximation may occur during inflation, leading to interesting 
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observational effects in the power spectra of cosmological perturbations. An archetypical 
model to study such violations is given by the potential 



1 arctan 

TT V mpi 



(6.51) 



introduced with = 5 in Ref. [6T]. In Fig. I6.13l we show that the MCE result expanded to 




-1.5 -1 -0.5 0.5 1 1.5 



Log{k/k^) 



Figure 6.13: Spectrum of the Mukhanov variable Q {Pq = |(5a;P/(2 vr^)) for the arctan 
model of Eq. fl6.5ip . The dots represent the numerical values, the solid line the analytic 
results from Eqs. fl6.39ap and fl6.39cp obtained by the second slow-roll order MCE approx- 
imation and the dashed line those given by the first order slow-roll approximation, k^, 
crosses the Hubble radius at (p^ ~ — 0.3mpi (the lines are normalized at lO'^^'^^k^) and the 
spectrum is evaluated at ^ 55 e-folds afterwards. 

second slow-roll order provides a very good fit for the power spectrum of scalar perturba- 
tions even in situations where the slow-roll parameters are not very small (see Fig. 16.141) . 
This example also shows that second slow-roll order results are much better than first 
order ones (analogous results were also obtained with the GFM in Ref. [68] ) . 
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Figure 6.14: Evolution of with the value of the inflaton (in units of mpi) in the arctan 
model: ei (solid line), €2 (long-dashed line) and €2/ H = e2e3 (short-dashed line). 
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Conclusions 



In this thesis, we have developed some approximations of the WKB type for the purpose 
of estimating the spectra of cosmological perturbations during inflation. In particular, in 
Section we found general formulae for the amplitudes, spectral indices and a- runnings 
of the fluctuations to next-to-leading order both in the adiabatic expansion of Ref. [37] 
and a new perturbative expansion which makes use of the Green's function technique. 
The most sophisticated Method of Comparison Equations was introduced in Section 14.41 
and then applied to cosmological perturbations in Section 15. 3[ By construction (i.e. by 
choosing a suitable comparison function), this approach leads to the exact solutions for 
inflationary models with constant horizon flow functions e^'s (e.g. power-law inflation). 
The main result is that, on using this approach to leading order, we were able to obtain 
the general expressions fl5.45ap - fl5.45cl) for the inflationary power spectra which are more 
accurate than those any other method in the literature can produce at the corresponding 
order. In fact, the MCE leads to the correct asymptotic behaviours (whereas the standard 
slow-roll approximation fails [67]) and solves the problems in computing the amplitudes 
which were encountered with both the standard and improved WKB method analised in 



In Section 16.11 we have applied our methods to power-law inflation in order to test it 
against exact results. It is known that the spectral indices and their runnings are obtained 
exactly at leading order in the WKB approximation, hence we focussed on the spectra 
to next-to-leading order. For the improved WKB method, we have found that the per- 
turbative corrections remain too small to yield any significant improvement, whereas the 
adiabatic expansion to next-to-leading order reproduces the exact amplitudes with great 
accuracy (see Ref. [62] for a similar result with the uniform approximation). This result 
does not however mean that the perturbative expansion will not lead to significant correc- 
tions in different inflationary scenarios. One way of understanding the difference between 
the two expansions in the power-law case may be the following. It is known that the 
Born approximation in quantum mechanics (analogous to our perturbative method with 
the Green's function) is good for high multipoles (angular momenta) ^ of the expansion 
in spherical harmonics. In our approach, the parameter (3 plays the role of the multipole 
index ^ for the energy levels of the hydrogen atom, as can be seen on comparing the 
frequency fl6.2l) with the expression given in Ref. [31] (see also Ref. [SB])- The approxima- 
tion with the Green's function is therefore expected to yield more significant corrections 
for large values of which coincide with the regime of fast-roll. Indeed, the leading 
order becomes more and more accurate for increasing \(3\. Since the interesting regime 



Sections 15.11 and 15.21 
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for inflation involves small values of /? ~ —2, it is instead the adiabatic approximation 
which seems better since the horizon flow functions evolve slowly and/or the states are 
quasi-classical B This may be the reason whereby the adiabatic approximation works bet- 
ter at the next-to-leading order. However, in more general cases, either or both methods 
may contribute significant corrections. In this respect, let us remark that the adiabatic 
expansion can be straightforwardly applied only for a linear turning point 0, whereas the 
perturbative Green's function method is not so restricted. 



In Section 16.21 we have shown that the improved WKB treatment of cosmological 
perturbations agrees with the standard slow-roll approximation [67] to within 0.1 %, finally 
resolving the issue of a 10 % error in the prediction of the amplitudes to lowest order which 
was raised in Ref. ^58j. The next issues are the inflationary predictions to second order 
in the slow-roll parameters. After the results on the running of the power spectra [TB] , 
second order results have been obtained using the Green's function method with the 
massless solution in a de Sitter space-time [75l EH] . Here we have employed the leading 
WKB approximation to obtain the scalar and tensor power spectra to second order in 
the slow-roll parameters. The key technical point is to use Eq. (16.211) . which allows one 
to express the power spectra on large scales as a function of the quantities evaluated 
at the Hubble crossing, but with no explicit dependence on AA^. As one of our main 
results, we find that the polynomial structure of the power spectra in the e^'s obtained 
from the WKB method agrees with the one arising from the Green's function approach 
of Refs. [751 [68] . In Section 16.2.31 we employ the next-to- leading WKB approximation 
to second order in the e^'s and, by requiring that the power spectra do not explicitly 
depend on AA^ (a property which instead was derived both in Sections 16.2.11 and 16.2.21) . 
and using the expressions found in Section 16.2.21 we obtain unique expressions for the 
spectral indices and runnings. Our findings show that different ways of approximating 
cosmological perturbations show up at the leading order in the amplitude, but in the 
next-to-leading order in the derivatives of the power-spectra with respect to wave number 
k. The accuracy in the theoretical predictions on spectral indices to second order in the 
slow-roll parameters, evaluated at A; = /c=k, is now striking: the Green's function method 
yields a coefficient C ~ —0.7296 which is replaced by -Dwkb* — —0.7302 in the improved 
WKB method, leading to a precision of 1 part on 1000 in the predictions of the coefficients 
of C(e^) terms in the spectral indexes. 

In Section Rj.2.3l starting from the general results fl5.45al) - fl5.45cp . we have also com- 
puted the full analytic expressions for the inflationary power spectra to second slow-roll 
order in Eqs. fl6.37ap and ( ]6.37bl) and found that the dependence on the horizon flow 
functions e^'s is in agreement with that obtained by different schemes of approximation, 
such as the GFM [75] and the improved WKB approximation (see Sections |6 . 1 1 and [6^ . 
Moreover, the results obtained with the MCE do not contain undetermined coefficients, in 
contrast with the second slow-roll order results obtained with the WKB*. Let us conclude 
by remarking that, just like the WKB approach, the MCE does not require any particular 

^Cosmological perturbations are amplified by inflation evolving from vacuum to highly squeezed states, 
which resemble classical states in the amplitude of fluctuations [50] . 
''The same limitation seems to affect the method used in Ref. [12] ■ 
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constraints on the functions e^'s and therefore has a wider range of apphcabihty than any 
method which assumes them to be small. As an example, we have discussed in some detail 
the accuracy of the MCE for the massive chaotic and arctan inflationary models. We have 
shown that the MCE leads to accurate predictions even for a model which violates the 
slow-roll approximation during inflation. 

In Appendix O we have also applied the Method of Comparison Equations to a dif- 
ferent problem, namely the scalar field dynamics on the Schwarzschild background. We 
have considered modes with energy lower than the peak of the effective potential in the 
corresponding wave equation, so that there are two zeros (turning points). The com- 
parison function has then been chosen of the Morse form with an argument which we 
have estimated analytically as well. This procedure was shown to be able to produced a 
fully analytical approximation of the wave-functions with good accuracy over the whole 
domain outside the Schwarzschild horizon. 
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Appendix A 

First Order Fluctuations 



We shall here summarise some technical aspects concerning the scalar, vector and tensor 
fluctuations of the geometry. As in the main text, we will assume a conformally fiat (i.e., 
K = 0) background metric of FRW type. 

First off, the fluctuations of the Christoffel connections can be written to first order 
in the amplitude of the metric fluctuations as 



2 

+\5g''' {dp ^'^gua + - ^'^g^p) , (A.l) 

where denotes quantities evaluated for the background geometry and the usual con- 
vention of summation over repeated indices is adopted. Analogously, the first-order fluc- 
tuation of the Ricci tensor is given by 

+s r^, (°)rf^ +(°) r-, sr^^-s rf^ ^r^, - ^ r^^ s r^, , (a.2) 

which leads to 

6R'; = Sg'-^^'^^R^.+^'^g'^'SR^, (A.3) 
Finally, the fluctuations of the Einstein tensor 

6g; = 6r;-^6';6r, (a.s) 

can be easily obtained from Eqs. flA.Sp and flA.4p and the perturbation of the covariant 
conservation equation is 

ST^"" +(°) r;;^ st^^ + s r;;^ (°)t°" +(°) r^^ 5 t"^ + 5 r^^ (°)t"^ = . (a.6) 

111 



112 



Appendix A. First Order Fluctuations 



The Christoffel connections of the flat FRW background are exphcitly given by 

while the components of the Ricci tensor and scalar are 

(°)i?oo = -3 7^', (X° = -4^' 

(o)i? = _A (n' + n') . (A.8) 
A.l Scalar modes 

The scalar fluctuations of the geometry can be written to first-order as 

2 

4 5-00 = 2 , 4 fl'°° = ^ 

6s g^, = 2 6i, - d, d, E) , 6, g'' = - E) 

Ss go^ = -a' d,B , s^gO^ = _]_0'B . (A.9) 

Inserting Eqs. (]A.9I) into Eq. ( lA.ll) . the first-order components of the (scalar) Christoffel 
fluctuations can now be easily computed 



Os J- 00 



4 r° = -[ip' + 2n{(j) + tfj)] + di d, {E' + 2nE-B) 

6,rio = d' {(p + B' + HB) 

4 T'^j = 8''^ 6ij -d^ij 5f - dj ^ Sf + di d^E-Hd'^B 

6sTi, = -^'6i + did' E' . (A. 10) 

The first-order fluctuations of the Ricci tensors are 

^3 i?oD = V [0 + (5 - E')' + n{B- E')] + 3 + 7^ (0' + ij')] 

s.Roi = di [{n' + 2n^)B + 2{^' + n^)] (a.ii) 
4 Rij = + 2{H' + 2 re) {^ + (j)) + n (0' + -v^t/j 

+nV^ {B-E')] 
+di dj [{E' -B)' + 2{n' + 2n^) E + 2n{E' -B)+ij-(f)] , 

while the first-order fluctuation of the Ricci scalar is 

s,R = ^ {34," + 6 (n' + n^) 4> + 3n {(p' + 3ij') 

[(j) - 2iP + {B - E')' + 3n{B- E')] } . (A.12) 



A. 2. Vector modes 
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Using Eqs. (lA.lip into Eqs. (1A.3P and (]A.4p the fluctuations of tlie Ricci tensor will be 



4i?o = l{v2 [(f) + {B-E')' + n{B-E')]+3[ip" + n{<i)' + ^') + 2n'<p]} 

(X 
(X 

-^d,d^ [{E'-B)' + 2n {E' -B)+ij-(P] 

s,Rl = ^d' [-{ij' + n^) + {n'-n') B] . (A.13) 

Finally the fluctuations of the Einstein tensor are 

= ^ {v^ij -nv^ {B - E') -sn ii^' + n4>)} (a.m) 
5,g{ = \ {-2i)" -2 [n^ + 2n') (p-2n<p' -AHi)' 

_ V2 [(p-^ + (B-E')' + 2n{B-E')]}6i 
-\did^ [{E' - B)' + 2n{E' - B) + ^ - ^] (A.15) 
4^° = 6R°i . (A. 16) 

A. 2 Vector modes 

The vector fluctuations of the geometry can be written as 

qi 

5v 9^J = -a' {d, E, + d, E,) , 6, g'^ = 1 {d' E^ + E') (A. 17) 

and the corresponding fluctuations of the Chrristoffel connections 

^vr°. = i {d,Sj + d„s,)+n {d,E, + d,E,) + ^ {d,E'^ + d,E'^ 

r^o = -s'' -ns' 
1 

2 



+ ^dj {d'E, + d,E'')+^d, [djE'^ + d'E,) 



rlo = -I {^^ - S.) + \ {d^ El + E^') (A.18) 
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where Eq. (I2.1ip for Fj and Si has been used. The fluctuations of the Ricci tensor are 

+n {d, f; + dj Fi) + ^ {di Fj + dj Fi) {2n' + A n^) , (a. 19) 
A. 3 Tensor modes 

Let us now considering the tensor modes of the geometry, that is the two polarisation of 
the graviton 

5t gij = -a^ hij , St g'^ = — . (A.20) 

a 



From Eq. flA.ip the tensor contribution to the fluctuations of the connections can be 
expressed as 

5t rj,. = i [d, hi + h) - K,] , (A.21) 



and the Ricci tensor becomes 

1 
2 



5, Rij = ]- [h'l^ +2nhij+2 [W + 2 n^) - V' hij] (A.22) 



(A.23) 



2a? 

A. 4 Energy- Momentum Tensor 

The background energy-momentum tensor of a perfect fluid is 

T^,u = {p + p) Uy -pg^,u ■ (A.24) 
From the normalisation condition 

g^,u''u'' = l, (A.25) 
one finds that Uq = a and 6u^ = — 0/a. Hence, the scalar perturbations become 



1 



5,T^ = 6p, 4T"'^ = - (5p-2p0) (A.26) 

a 



A. 4. Energy-Momentum Tensor 
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and 

6,T^ = -6p6i (A.27) 

4 = ^ [6p6'^ + 2p {tp 5'^ -d'd^ E)] (A. 28) 

5s = {p + p)v' (A.29) 

4T°' = ^[{p + p) v' + d' B] . (A.30) 

(X 

where we have defined Su^ = v^/a. The velocity field can be split into divergenceless and 
divergencefuU parts 

v' = d'v + V, d^V' = (A.31) 

For collisionless matter, the anisotropic stress must be considered both in the Einstein 
equations and in the covariant conservation equations. The anisotropic stress is then 
introduced as 

6sT^ = -6p6i + I^ , (A.32) 
where 11^ = T/ — 61 For the scalar fiuctuations we adopt the notation 

dj d' ni = - (p + p) VV . (A.33) 
The perturbed velocity field can be defined in two ways, 

4 n = {p + p)uo6u' = {p + p) v'^{p + p) d' V (A.34) 

or 

5s Tf = {p + p) uHui = {p + p) v, = {p + p)div , (A.35) 

where, from the normalization condition flA.25p . Uq = a and vP = 1/a. The velocity fields 
V and V defined in Eqs. ( 1A.34I) and ( ]A.35I) are not equivalent. In fact 5ui = avi and 
— yi- jo^^ Using 

bu' = 6 {g'^'u^) =6g'°uo + g'''Suk (A.36) 
S g^o = -d'B/a^ , (A.37) 

and inserting the definitions of 6u^ and 6uk in terms of f * and Vk we have 

Vi = -v,-d,B . (A.38) 

The gauge transformations for Vi and Vi are correspondingly different, 

Vi ^ Vi = Vi + di e' 

Vi^Vi = Vi - di eo . (A.39) 
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Appendix A. First Order Fluctuations 



In this thesis we use the velocity field as defined in Eq. (lA.34p . 

Let us now consider the energy-momentum tensor of a scalar field (p characterized by 
a potential V{(p), 



Denoting with x the first-order fiuctuation of the scalar field ip, we have 



-g^^d^ipdpip-V 



(A.40) 



(A.41) 



- Ss da<pdp<p + g"^ daXdp<p - j^X 



or, exphcitly, 



X V + 2 y + a' 



dV 

dip 



X 



ip'diX + di B 
6. 



2a2 



V 



xa" + ilj) ip'^ + 2a'Vij 



+2 a' 



/2 



2a2 



dip 



- V 



di dj E . 



(A.42) 



On using 



00 



/2 



(A.43) 



the perturbed components of the energy-momentum tensor can then be written as 



that is 



"S -'-I 



-<py^'' + ^a^X + x'^' 
dip 

,2 dV ^ , A 

+^<^X-X^]Oi 
dp J 



±^'d^X-^d^B . 



(A.44) 

(A.45) 
(A.46) 
(A.47) 



The covariant conservation of the energy-momentum tensor implies the validity of the 
Klein-Gordon equation which can be written as 



dip 



dV 



(A.48) 



A.5. Gauge invariant scalar perturbations 
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From Eq. (1A.48P the perturbed Klein-Gordon equation is 

5s [d^dpip-^'^^Vlpd^if] 



(A.49) 



or, more explicitly, 



+20^a2-<^' (</)' + 3 ^') + <^'V' = . 



(A.50) 



A. 5 Gauge invariant scalar perturbations 

We show here the generalised evolution equations for scalar fluctuations obtained without 
fixing a particular gauge. From Eqs. (lA.14p . (lA.26p . (1A.16P and (1A.29P the Hamiltonian 
and momentum constraints become 



V^ip-nV^ {B-E')-3n i-^' + n<f)) = AnGa'^6p 
d' [{^' + n(j)) + {n^-H') B] = -AnGa^ {p + p) v' 



(A.51) 
(A.52) 



From Eqs. (lA.lSp and (1A.28P we obtain the spatial components of the perturbed equations. 



atxGo' [5p5i - n^] = [ip" + {n^ + 2n') (P + n {<p' + 2 v^')] 

+ \ V [0 - ^ + (5 - E')' + 2n{B-E' 
-\ di [<p-i) + {B- E')' + 2n{B-E' 



(A.53) 



From Eq. flA.SSp we can separate the traceless part as 



AirGa^dp = ip" + H {(f)' + 2'iIj') + {U'^ + 2n') (j) 

+ lv2 [^-^ + ^B-E')' + 2n {B-E' 

o 

UtxGo^ {p + p) a = [(l)-ip + {B - E')' + 2n (B - E')] . 
The (0) and (z) components of Eq. flA.6p can be written 
= 9o 5s T°° + d, 5s + (2 5s r°o + 5s T^o) 

+ (2 (°)r°o +(°) F^o) 5s T°o +(0) r°. 5s + 4 r° wr^- 

= 9o 5s TO^- + 9fc 5s T'^^ + (5s + 5s F™,) (o)T'=^- 

+ ((°)F°o +(°) F^o) 5s T°^- + 5s F^oo + 5s F^^ (°)T^- + 2 (°)F^,, 5s T 



(A.54) 



(A.55) 
(A.56) 



OA: 
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Appendix A. First Order Fluctuations 



Inserting the perturbed connections (lA.lOp into Eqs. (1A.55P and (1A.56P we have 



= S p' - 3ij' {p + p) + {p + p) e + 3n {Sp + 6p) + {p + p) V^E' 
^ = {p + p) e' + d [p' + p' + AH {p + p)] + {p + p) V' B' 
+ [p' + n {p + p)]V^ B + V^5p+{p + p)V^<P , 



(A.57) 
(A.58) 



where the divergence of the velocity field 6 = diV"^ = did^ v has been introduced. 



Appendix B 

Slow-Roll vs Horizon Flow 



In this Appendix, we compare the HFF we used in the text with some other (equivalent) 
hierarchies. 

In Table IB. 11 we give the relations between the slow-roll parameters defined by some 
other authors and the HFF defined in Eq. fl3.65p . 



J. Martin et al. [69j"f 


E.D. Stewart et al. [67, 75j* 


S. Habib et al. [73] 


S. Habib et al. [621 EH 


e = ei 


e = ei 


5 = ei-\t2 


(^1 = i £2 - ei 


^2 = 2 ^2 £3 


(52 = -fei e2 + 2e? + iei + ie2e3 



Table B.l: Relation between (some) slow- roll parameters and the HFF used here and in 
Ref. [58j. *In Ref. [67j ei, 5, and (p 5/H cj) are used instead of e, 5i and 82 respectively. '''In 
Ref. |69j ^ is used instead of ^2- 



We also compare the Hubble-slow-roll (HSR) and the potential-slow-roll (PSR) pa- 
rameters (see Ref. [HD [E2j) with the HFF. We recall that some authors redefine and 

as and and write only e, rj and i^, which can be confusing. The relations between 
the ej's and HSR parameters are 

1 

= ei , Vh = - -^2 

(B.l) 

^1 = - o ei £2 + - £2 63 , 
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Appendix B. Slow-Roll vs Horizon Flow 



and for the PSR parameters we have 

2 



1 + 



62 



2(3-ei: 



6 61 - I 62 - 2 + I ei 62 - i 62 - I €2 £3 

3-61 



261 --62 



(2 61 - 6 - 62) [8 6? -6 6? (4 + 3 62) + 61 62 ( 18 + 6 62 + 7 63) - 62 63 (3 + 62 + 63 + 64) 



4 (3 -61)^ 



~4 6i - 3 6162 + -6263 



(B.2) 



where we have shown both the exact formulas and approximate relations to the order of 
interest . 

We complete this review with the first three exact connection formulae between the 
PSR and HSR parameters (EH ES] 



ev (3 - iihY 



(3 - en 



,2 ' 



Vv 



3{eH + VH)-vl-eH 

3-6// 



(B.3) 



^ 9 (36//r/// + elf - envl) ' 3 (4////^!/ - vUh + ^h) + Vh(t% 

(3 - 6//)2 



and so on with an the next HSR parameter. 



Appendix C 

An application of MCE 



Here we use the Method of Comparison Equations to solve the differential equation (14.321) 
with the frequency uj^{x) = — 1 + a sin(x), which does not admit an exact solution. In 
fact we shall solve in a semi-analytical way the equation 

= -^[e^-l + a sin(x)] x{x) . (C.l) 

We shall give the semi-analytical results to leading and next-to-leading order of our ap- 
proximation, correspondingly to O*'^ and l^*-order in e. To build the semi-analytical solu- 
tion we need the first two terms of the asymptotic expansion (14.381) . i.e. 

yix) = Voix) +eyi{x) , (C.2) 

and also the function 

The function y{x) appears in Eq. (I4.34p as a prefactor (See Eq. (14.361) ). while (t(x) emerges 
in the argument of the analytical part of the solution (14.341) . i.e. in U [a{x)]. The above 
quatities will be estimated numerically taking e = 1 and a = 1/100. The initial conditions 
are u{—20) = 1 and «(— 20) = 0, while the range of the plots is —20 < x < 4. We recall 
that the quantity o"o changes from first to the second order of approximation, as we have 
shown in section l4.4.1l for the linear case. We shall use do = for the leading order, while 
at the next-to-leading order we shall recalculate its value utilising the previous order. 
Finally, the analytical part will be a general solution of Eq. (I4.33p . 

In the following we shall consider two comparison functions; a linear one, i.e. 0^(0") = 
cr, and a more complicated one, i.e. 6^(cr) = e*^ — 1 



C.0.1 Linear comparison function 

For the hnear (i.e. Langer) case yo{x) and yi{x) are given by Eqs. (I4.5ip and (I4.55ap . 
respectively. In Fig. lC.ll we show the function yi{x), while in Fig. lC.2l we plot the function 

^We use a "small" value for a in order to avoid other turning points in the oscillating region. 
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Appendix C. An application of MCE 



y{x) at 0* and l^*-order in e. The function cr{x) at 0* and I'^^-order in e, calculated 




Figure C.l: Function yi{x) for a linear comparison equation. 




-20 -15 -10 



Figure C.2: Function yo{x) (solid line) and y{x) = yo{x) + e yi{x) (dashed line) for a linear 
comparison equation. 

from Eq. fIC.Sp , is shown in Fig. IC.3I 

The analytical part of the solution follows from 



and is a linear combination of Airy functions Ai and Bi 

U [a{x)] = A-Ai [e-^/V(x)] +B-Bi [e-^/^a{x)] . 
Finally the complete solution becomes 

X{x) = y{x) {A ■ Ai [e"^/V(x)] + B ■ Bi [e'^/^a{x)] } 



(C.4) 



(C.5) 



(C.6) 



In Fig. lC.4l we exhibit the exact numerical solution and the first and second order solutions 
in our approximation scheme. 
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Figure C.3: Function a{x), for a linear comparison equation, to O*'^ and I'^^-order in e 
(solid and dashed line respectively). One can not distinguish the two curves. 




Figure C.4: Solution of Eq. ( IC.ll) for a linear comparison equation; the solid line is 
the full numerical solution, while the dashed and dotted ones are the O*'^ and I'^^-order in 
e respectively. 



CO. 2 A more complicated comparison function 

For the comparison function O^(cr) = e'^ — 1 the quantities yo{x) is implicitly given by 
Eq. (HZHD, while yi{x) is given by Eq. (HT9ll on using Eqs. f!4.8()al) and (HTTHIl . In Fig. [031 
we show the function yi{x) while in Fig. IC.6l we plot the function y{x) at 0**^ and I'^^-orders 
in e. In Fig. IC.7I we exhibit the function G{x) from Eq. fl4.78l) . The function a{x) at O*'^ 
and I'^^-orders in e, calculated from Eq. (IC.Sp . is shown in Fig. IC.8[ 
The analytical part of the solution follows from 

and it is a linear combination of Bessel functions / 
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-0.005 



Figure C.5: The function yi{x) for a more complicated comparison function. 




. 97 



-20 -15 -10 



Figure C.6: The functions yo{x) (sohd hue) and y{x) — yo{x) + eyi{x) (dashed hne) for 
a more comphcated comparison function. 



5 X 10" 
4 X 10^ 
3 X 10"^ 
2 X 10*^ 
1 X 10^ 
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Figure C.7: Function G{x). 



Finally the complete solution becomes 



X{x) = y{x) { A ■ I_2i 



2^ 



(C.9) 
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Figure C.8: The function a{x), for a more complicated comparison function, to 0**^ and 
l^t.orders in e, solid and dashed lines respectively. As in the linear case one can not 
distinguish the two curves. 



In Fig. IC.9I we exhibit the exact numerical solution, the first and second order solutions 
in our approximation scheme. 




X 



Figure C.9: Solution x{x) of Eq. (IC.ip for a more complicated comparison function; the 
solid line is the full numerical solution, while the dashed and dotted lines are respectively 
the O*'^ and l*^*-orders in e. 



In Fig. IC.4I we can note that the leading solution calculated with a linear comparison 
function is close to the true solution in the oscillating region, while the next-to-leading 
solution is better in the exponential region. It is quite natural because an iterative scheme 
is more stable in regions where a function under consideration has a monotonic behaviour. 
From Fig. IC.4I and Fig. IC.9I we can see that the more complicated comparison function 
leads to better solutions than the linear one. This feature is true because the more 
complicated comparison function remains close to the original frequency in the whole 
range of interest, in contrast with the linear case. 
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Appendix D 

Other results in the literature 



We use Table IBTTI and Eqs. (IB.ll) . to compare the results for spectral indices - and for PS 
and runnings when available - in terms of the HFF in Tables ID.ll and ID.2I 

In Table ID.ll we summarize the results for scalar and tensor perturbations given in 
Ref. [68j, as calculated with the method proposed in Ref. [75]. We also refer to some other 
papers which give the same results, although expressed in different hierarchies. 

Let us then compare with the results obtained by the uniform approximation [73l l62l 
The integrand in Eq. (9) of Ref. [73j can be obtained from our Eq. (15. ip by the 
change of variables a; = In {—ki]), u = e~^/^ /i, which of course differs from our Eq. fl5.9l) . 
Therefore, the integrands in our Eq. ( 14.14b[) and Eq. (9) of Ref. [73] differ by 0{ef) terms, 
once both are expressed in conformal time. We also note that in Refs. [73t l62| [71], the 
authors just give expressions for the spectral indices which depend on k (i.e. without 
referring to a pivot scale k = k^). Further, it is not easy to see that they obtain the 
correct scaling in In (k/k^,), which yields the a-runnings, and the ej's are evaluated at the 
classical turning points for the frequencies. 

In Tables ID. 21 we display the results for scalar and tensor perturbations given in 
Ref. [H]. As it was properly stated in Ref. [71], the results for scalar and tensor spectral 
indices given in Ref. [731 ES] were not fully expanded to second order in the slow-roll 
parameters. 
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Appendix D. Other results in the literature 



Table D.l: Results obtained with the Green's function method [68] suggested in Ref. [75] : 
C = ln2 + 7£; — 2 ~ —0.7296 with 7^; the Euler-Mascheroni constant. We display the 
original results for spectral indices, a-runnings and tensor-to-scalar ratio evaluated at the 
pivot scale k^,. 



E.D. Stewart et al. [6^, E.D. Stewart et al. [75] , 
S.M. Leach et al. [68], A.R. Liddle et al. [811182]. 



Vc 



TV €1 m-p 



|l-2(C + l) ei-Ce2+ (2^2 + 2^ + ^-5) e? 



+ [-2ei-e2 + 2{2C + 1) el + {2C-l) €162 + C ej - C 6263] In (^ 



+ 1 (4e? + 2eie2 + ei-e2e3) In^ (^) } 
n = ^ {l-2(C + l) ei+ (2^2 + 2^ + ^-5) e2+(-C2-2C+g-2) 6162 
+ [-2ei + 2(2C+l) el-2{C + 1)6,62] In (a) + 1 (4 _ 2 ei £2) In^ (a)| 

ns-l = ~26i- 62 - 26i-{2C + 3) 6162-C6263 , nx = -2 d - 2 e? - 2 + 1) d €2 
as = -2 ei €2 - 62 63 , ot = -2 ei €2 



/?= 16ei 



1 + Ce2 + (C-^ + 5) eie2 + (|C2-^ + l e2+ ic'S.i., ,^^3 



Table D.2: Results obtained with the uniform approximation. 

S. Habib et al. [74j 
ns-l = -2ei-e2-2e2 - (f -27r) 6162 -{f-n) 6263 
nT = -26i-26j -(T-f^) ^1^2 



Appendix E 



Undetermined parameters in the 
WKB* 

To clarify the impact of the undetermined parameters bs, &t, and c/s = 2 on the results to 
second order in the HFF, the percentage differences between the numerical coefficients of 
some second order terms in the power spectra are displayed in Table lE.ll and of some of 
those in the scalar-to-tensor ratio in Table [K2i All the entries are for 63 = = c^s = 2, as 
is used in the main text, since other choices would give larger differences, as we explicitly 
show in Fig. lE.ll for the numerical coefficient in front of ei €2 in Ph (since the behaviour 
is the same for all coefficients, we do not display other plots). 




1 1.5 2 2.5 3 



Figure E.l: Percentage difference between the numerical coefficients of 6162 in from 
WKB* and GFM as a function of 6t- The difference vanishes for 6x — 1-93 2, the 
value used in the main text). 
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^I'l Ph 


4;P, 


eie2;n(^s = 2) 


ei £2; Ph{bT = 2) 


WKB 


2.3% 


3.0% 


353.3% 


37.2% 


WKB* 


0.03% 


0.06% 


29.1% 


2.0% 



Tabic E.l: Percentage differences between the numerical coefficients of some second order 
terms in the power spectra given by WKB and WKB* with respect to those obtained by 
the GFM. 





£263; {ds = 2) 


ei £2; {bs - 6t = 0) 


4 


WKB 


56.7% 


1.5% 


117.3% 


WKB* 


1.3% 


0.03% 


1.5% 



Table E.2: Percentage differences between the numerical coefficients of some second or- 
der terms in the scalar-to-tensor ratio given by WKB and WKB* with respect to those 
obtained from GFM. Note that the numerical coefficient of €2 €3 in R is the same as that 
in Pq. Further, the coefficient of ei €2 does not depend on any undetermined parameters 
if 6s = 



Appendix F 

Perturbative approximations 



In this Appendix we present a possible method of obtaining the complete result to second 
order in the HFF. We start by noting that any given function / = f{x) = f{e,x) can be 
written as the power series 



n=0 



X 

n\ 



f{ei,x) = V ^/„(x) 



n=0 



(F.l) 



where e, = e the HFF evaluated at the horizon crossing (i.e. at a; = x*) and 

/(ej,x) = /(ej = ei,x). On using Eqs. fIF.ip to express both u"^ and x i^i Eq. (14.61) . we 
obtain 



n=l 



H2 °° 



■m=0 



n=l 



(F.2) 



where we have used 



dx^ 



+ u}^{x) 



Xo{x) = 



(F.3) 



with co'q(x) = uj'^{ei = ej;x) and Xo{x) is a linear combination (i.e. the MCE leading 
order solution) of Eq. (15.421) for B^(x) = uJq{x) and cr — >■ x in the argument of the Bessel 
functions. One may solve Eq. (IF. 21) for higher orders by introducing 



nl 



(F.4) 



and using the standard formulae for second order differential equations with a source 
term, so that 



gn[x) = - 



xliy) 



n—l 



m=l 



n — m 



(F.5) 



The problem has therefore been reduced to computing 2 n integrals. 



131 



132 



Appendix F. Perturbative approximations 



The above formalism can be modified in order to provide an alternative way for cal- 
culating the power spectra, spectral indices and runnings to second order in the HFF. We 
start from Eq. (15.11) using the exact relations 



3 1 1 2 1 

ei + - £2 - - ei £2 + - €2 + - €2 eg 



With a redefinitions of the wave-function and variable 

ijj = e~^^'^ /i , x = In {—k 7]) 



we obtain 



where 



dx^ 



^{x) = 



(F.6) 



(F.7) 
(F.8) 



~,2 



,2£ 



,2£ 



r]'a'H' (2 - d) 



^ - r]^a^H^ f 2 - ei + ^ 62 - ^ ei e2 + ^ + ^ e2 eg 



On now using the approximate relation 

1 



r] 



aH 



(l + ei + e2 + eie2) 



the frequencies become 



ul ~ e^* 



- + 3ei +4ef + 4 6162 



9 o 3 , 2 13 1 2 1 

- + 3 6i + - 62 +4 6i + y 6162 + 4 ^2 + - 6263 



(F.9) 



(F.IO) 



(F.ll) 



where we only keep the second order terms in the (still x-dependent) HFF. We can now 
expand linearly each HFF around the horizon crossing to obtain 



d^T 
dx^ 



d^ 
dx^ 



e2-_(^ + 3ei + 4e? + 4eie2 



Tp'j: = —3 ei €2 xipT 



(F.12) 



^2x 



9 . 3 „ .2 13 _ 1 . 1 _ 
- + 3ei + -e2 + 4ei + y £162 + 4 ^2 + 2 ^2^3 



-3 1 ei £2 + - 62 £3 ) X ?/'s , 
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which, beside the use of x instead of x and the appearance of uj instead of cj, is equivalent 
to the expansion in Eq. (IF.ip for / = cD^ and can be solved by means of a standard 
perturbative approach (such as in Ref. [75]). The method is close to GFM and gives the 
same results of Eq. (41) in Ref. [75] and Eq. (29) in Ref. [68] for the scalar and tensor 
power spectra, respectively. 

Let us finally note that, for the GFM, the scalar spectral index and its running are 
calculated by taking the derivative of Eq. (41) in Ref. [7S], i.e. the power spectra evaluated 
at A; = a if, in terms of In k, while in our MCE (but also in the WKB approximation [701 
171] ) the dependences in ln(A;/A;=K) are obtained explicitly (see, for example, Eqs. f l6.37al) 
and (lOTbl) ). 
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Appendix G 



MCE and 



Black Holes 



In this Appendix we employ the Method of Comparison Equations to study the propa- 
gation of a massless minimally coupled scalar field on the Schwarzschild background as 
an example of a potential with two turning point in the wave equation. In particular, we 
show that this method allows us to obtain explicit approximate expressions for the radial 
modes with energy below the peak of the effective potential which are fairly accurate 
over the whole region outside the horizon. This case can be of particular interesting, for 
example, for the problem of black hole evaporation. 

Wave equations on black hole backgrounds cannot be solved exactly in terms of known 
simple functions (for a review, see, e.g., Ref. [EH]). Although it is usually fairly easy 
to find asymptotic expressions near the (outer) horizon or at spatial infinity, numerical 
analysis is what one eventually employs to obtain a detailed profile of the wave-functions. 
However, our understanding of the field dynamics on curved backgrounds would still 
benefit of analytical approximations to the exact modes, for example in the study of the 
(renormalized) energy-momentum tensor of quantum fields [33] and, in particular, of the 
Hawking radiation [Ml [85] , or of quasi- normal modes [86] . 

By making use of the black hole symmetry, one can usually separate the wave- functions 
and reduce the problem to the one task of solving a second order differential equation for 
the radial part. The latter takes the form of the one-dimensional Schrodinger equation 
for the transmission through a potential barrier [S3] 



and, for the general case of a Kerr-Newman black hole, is known as the Teukolsky master 
equation for waves of spin- weight s = 0, 1 and 2 [ETJ [83]. 

In the next Section, we shall briefly review the equation of motion for massless scalar 
fields on a spherically symmetric black hole background. This includes the simplest case 
of linear perturbations of a Schwarzschild black hole which will be used to show the 
application of our method in all details (including the use of a convenient radial coordinate 
and rescaling of the wave-function). In Section [0.21 we shall then apply the MCE and 
obtain approximate expressions for scalar wave-modes with energy below the peak of the 



d^ 



+ Qs{x) V = , 



(G.l) 



dx^ 
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potential. This case is of particular interest, for example, for evaluating the grey-body 
factors involved in the Hawking effect [Ml [S3] . It also presents a major complication with 
respect to the typical case of Cosmological Pertubations, namely the presence of two zeros 
of the potential ("turning points"). We shall use units with c = G = 1. The last three 
Sections will give some tecnical results which complete this Appendix. 



G.l Scalar field on Schwarzschild background 



Let us begin by recalling that the metric in a static and spherically symmetric vacuum 
space-time can be written in general as 



ds^ = h{r) 



dr^ 
h{r) 



(G.2) 



where dQ is the area element on the unit two-sphere. The propagation of massless, 
minimally coupled scalar particles in this background is governed by the Klein-Gordon 
equation ^'^.^ = 0. Due to the symmetry of the background, the field $ can be de- 
composed into eigenmodes of normal frequency u and angular momentum numbers i, m 

as [83J 



$(t,r,fi) =e-*'^*r^"(n)i?(r 



(G.3) 



where are spherical harmonics. The Klein-Gordon equation then separates and the 
dynamics is described by the function R which satisfies the radial equation 



d 
dr 



h{r] 



dR{r) 
dr 



+ 



u"^ r^ 



_h{r) 



R{r) = 



(G.4) 



where i{£ + 1) is the separation constant. No exact solution of the above equation is 
known even for as simple a background as the four- dimensional Schwarzschild space-time 
with 



h{r) 



(G.5) 



where we have introduced the dimensionless radial coordinate f = r/rn and rn is the 
Schwarzschild radius. 

It is now convenient to introduce the "tortoise-like" coordinate 



dx 



dr 
r h{r) 



and the new radial function 

R [r{x)] = exp 



h{x') d x' 



(G.6a) 



(G.6b) 



similar coordinate was already used in Ref. 
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so that the radial equation takes the Schrodinger form 



1 



+ <^ Cu' f\x) + h[f{x)] - - h'[f{x)] - - 



1 d/;,[r(a;)] 



dx 



X{x) = . (G.7) 



where we also defined the dimensionless energy uj = coTh- For the four- dimensional 
Schwarzschild case we have 



x{r) = ln(r — 1) 



(G.8) 



so that the new variable x ranges from x(fH = 1) = ~cxo to x{r = +oo) = +cxd and x and 
x{x) are of Langer's form |37j. The Klein-Gordon equation finally becomes 



dx" 



+ uj^{x) 



X{x) = 



(G.9) 



where the "frequency'' 



^,2 iii+l)e^ e^(2 + e^) 



1 -Fe^ 



A (l + e" 



\2 ' 



(G.IO) 



is not necessarily a positive quantity. In fact, the above expression shows two qualitatively 
different behaviors depending on the values of uj and i. In particular, for a given angular 
momentum there exists a critical value 



{G.ll] 



such that if the energy u > 0)^ there are no turning points, otherwise there are two, say 
Lj{xi) = uj{x2) = with Xi < X2- The latter case is the one we want to analyze in detail 
(see Sec. IG.2.1l for more details and the solid line in Fig. IG.ll for an example). 



G.2 MCE solutions 

Referring to Section 14.41 for the basic formulae characterising the MCE, we can find the 
approximate solution of Eq. flG.9p . 

Xmce(x) = (^^^ ^ Uia) , (G.12) 

which is valid in the whole range of x and including the turning points. 

We are now dealing with a problem of the form (1G.9P with two turning points, Xi < 
X2- In order to implement the MCE, we need a "comparison frequency" with the same 
behavior. We shall use the Morse potential [89] 



e^a) = Ae^" - Be" + D 



(G.13) 
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The coefficients A and B are then fixed by imposing that the turning points of the 
comparison function are the same as those of the exact frequency, that is 



e((Ti = Xi) = 0((T2 = X2)=0 

and the coefficient D will then follow form the relation 





J Xl J X\ 



The system of non-linear equations (IG.14ap and (IG.14bp thus yields 



A 



{wi + W2 - 2 ^Wi W2) ' 



(G.14a) 
(G.14b) 

(G.15a) 



B = {wi+W2)A (G.15b) 
D = W1W2A , (G.15c) 

where the exact expression for ^ is given in Sec. IG.2.21 A sample plot of both the exact 
frequency and the comparison function thus obtained is given in Fig. IG.ll 




Figure G.l: Comparison between the exact potential —00"^ (x) (solid line) and the approxi- 
mate Morse potential — 0^[cr(x)] obtained after solving Eqs. (lG.14al) and (lG.14bl) (dotted 
line) for uj = 1/2 and i = 2. The two turning points are xi ~ —3.13 and X2 ~ 1.23. 

The "comparison solution" is a linear combination of confiuent hypergeometric func- 
tions [M] of the type iFi (for more details, see Ref. [SU]). 



Uia) 



D 



B 



l + 2iVD, -2iVAe' 



B 



2VA 



2^fA 



l-2iVD,-2iVAe'' 



(G.16) 



and o"(a;) is given implicitly by Eq. (14.311) with 

v/e%^dp = vW) 



(G.17) 



In 



2Ae'' - B + 2^/Ae^^{a] 



B 
2 VA 



2De-'' - B + 26-" V7J02(^ 
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The above expression is rather complex. Hence, we just consider the asymptotic relation 
between a and x to the left of the smaller turning point (that is, for x <^ Xi) 



a{x) ^ 

and to the right of the larger one (x ^ X2) 

B 



D 



X 



(G.18a) 



2Vv4 



(Tlx) ~ C.J e + x) . 



(G.18b) 



As an approximate expression for 0"(x) for x < /ti Xi (with /ti > 1) we shall then 
use Eq. ( IG.lSal) and the condition that o"(xi) = xi. Analogously, for x > ^2X2 (with 

^ 1) we shall use the approximate solution of Eq. fIG.lSbp to next-to-leading order for 
large x. Moreover, in the region between /ti Xi and ^2 X2 we shall use an interpolating 
cubic function. The reason for introducing two new parameters Ki and K2 ^ 1 is that the 
asymptotic forms in Eqs. (IG.lSal) and (IG.lSbl) may differ significantly from the correct 
ct(x) around the turning points Xi and X2 and suitable values of ni and K2 usually improve 
the final result. To summarize, we have 



X + Xi I 1 



Xq < X < KiXi 



Co + Ci X + C2 X^ + C-s X^, KiXi < X < K2X2 



(G.19) 



X + In 



X > K2X2 



with Ki and ^2 > 1 and Xq < Ki Xi is the value of x at which we wish to impose the initial 
conditions for x{^) ^"^^ its derivative @. The explicit expressions for the coefficients Co, 
Ci, C2 and C3 are rather cumbersome and will be given in Sec. IG.2.3[ 

In the upper left panel of Fig. IG.21 we plot cr(x) for the same values of u and £ used 
in Fig. IG.ll The solid line represents the numerical solution of Eq. f l4.3ip and the dashed 
line the approximate expression ( IG.IQI) with /ti = 1 and H2 = 3. It is clear that these 
values of Ki and K2 already lead to a very good approximation for a around the turning 
points Xi and X2 and that larger values of Hi or k,2 would not improve significantly the 
final result. In the upper right panel of Fig. IG.21 we compare an MCE approximate 
wave-function with the exact solution (solid line) obtained numerically for the same case. 
The MCE solution shown by a dashed line is given by the expression (1G.12I) with U{a) 
as in Eq. ( 1G.16I) and cr(x) determined by solving Eq. fl4.3ip numerically. The coefficients 
C± are fixed by imposing that the MCE mode and its derivative equal the corresponding 
values of the chosen numerical solution at x = —20. Since the two curves coincide almost 



^Of course, if we are interested in setting the initial conditions at xq 3> X2 the expressions must be 
adjusted correspondingly. 
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Figure G.2: Upper left graph: The function a = a{x) obtained by integrating Eq. (14.310 
numerically (solid line) and its approximate expression (1G.19P with ki = 1 and K2 = 3 
(dashed line). Upper right graph: Comparison between the exact solution (solid line) and 
the approximate MCE solution (1G.12I) with cr(x) evaluated numerically (dashed line) for 
C+ ~ 0.517 — 0.228 i and C_ ~ 0.517 + 0.228 z. Panel (a) shows the percentage error 
for the same approximate function a given in the upper left graph. Panel (b) shows the 
percentage error for the approximate x the upper right graph for x < Xi and panel (c) 
for X2 < X. All plots are for = 1/2, £ = 2 as in Fig. IG.li 



everywhere, it is clear that solving the approximate equation (I4.3ip for (j{x) is, to all 
practical extents, equivalent to solving the original equation (lG.9p . 

In order to further clarify this point, in panel (a) of Fig. IG.21 we plot the relative 
difference (in percent) between the approximate expression ( IG.IQP and the exact numerical 
solution of Eq. f l4.3ip . We can see that the error is less than 10% almost everywhere, 
except near the zero of a{x). In the same figure, we also show the percentage error for 
the approximate solution shown in Fig. IG.2I on the left of Xi and on the right of X2- The 
latter error remains small, too, except around the zeros of x due to the slight difference 
in phases. 

Finally, in Fig. IG.31 the dashed lines represent the same MCE solution with a{x) given 
by its analytical approximation in Eq. ( IG.IQP with ki = 1 and = 3. This is the fully 
analytical expression obtained from the application of the MCE and is still remarkably 
good, with a more significant error around the larger turning point {x2 ~ 1.23). From the 
three graphs in Fig. IG.3I and the upper left panel in Fig. IG.21 it clearly appears that such 
an error develops in the region between the two turning points, where the discrepancy 
between the exact a{x) and the approximate expression (1G.19P is indeed larger. In this 
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Figure G.3: Comparison between the exact solution (solid line) and the approximate 
MCE solution (1G.12P with a{x) given in Eq. (1G.19P with ki = 1 and K2 = 3 (dashed line) 
for u! = 1/2 and i = 2 (same as in Figs. IG.ll and IG.2p . Panel (a) shows the result for 
X < KiXi, panel (b) for kiXi < x < k,2 ^2 and panel (c) for x > ^2X2- 



respect, let us note that, had we imposed initial conditions at x ^ X2, the larger error 
would have occurred around the smaller turning point (xi ~ —3.13). Moreover, when we 
used a linear (or quadratic) interpolating function for niXi < x < 1^2X2, the accuracy 
was (expectedly) worse. Should one need better accuracy, a higher order polynomial 
interpolation between ni xi and K2 X2 could therefore be used instead. 

Let us now make a few remarks. Firstly, the application of the MCE is not as straight- 
forward as the standard WKB approximation. In fact, the MCE usually requires solving 
some technical tasks specific to the problem at hand in order to obtain the argument 
a = a{x) of the known function U = U{cr). In the present case, with two turning points, 
all expressions become more involved. The best approximation for cr(x) we found conve- 
nient to display in Eq. (1G.19P makes use of a cubic interpolation over the region containing 
the two turning points and is given in terms of two arbitrary parameters Hi and k,2 which 
can be fixed in order to minimize the error. Of course, one could employ higher order 
interpolating functions and further improve the result at the price of complexity if bet- 
ter accuracy is needed. Our aim for the work presented in this Appendix was mainly 
to test the effectiveness of the MCE to produce approximate wave modes on black hole 
backgrounds. 

In the next three technical Sections we give explicit expressions for the turning points, 
the calculation of the quantity ^ and some details of the cubic interpolation used to obtain 
a{x). 



G.2.1 Turning Points 

If we define w = e^, we see that C(j^(x) = for the exact frequency in Eq. (IG.lOp becomes 
the quartic polynomial equation 



+ 4:W^ + bw^ + cw + 1 = , 



(G.20a) 
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with 



6 
4 



1 

2^ 



Its four roots can be written as 



Wi 



where i = 1, 



-l + -W4--6 + ^ 
2 V 3 3 



2 



46- 16 -2c 



7 



4 and 



a 
7 



The coefficients ei = £2 



12 + 6^-120 

2 36 6c + 270^ + 432 - 72 6 



23 a 



/5+ v//?2-4«= 



-£4 



^1 



2i 



-00 = 0a 



(G.20b) 
(G.20c) 



(G.21a) 



-04 



(G.21b) 
(G.21c) 

(G.21d) 
1, and the values of 



Cj and £ determine the signs of the roots. There can be at most two positive roots, in 
which case, say, W3 < W4 < < wi = e^^ < W2 = e^^, and only xi and X2 can be valid 
turning points. The critical value of Cj in Eq. ( IG.lljl is then determined by the condition 

Q <Wi = W2- 



G.2.2 Evaluation of ^ 



In order to evaluate the integral in Eq. ( lG.14bl) . we start from the general relation (again 
with w = e^) 



^ /•ln(w) 



l {w -W3){W -Wi){w -W2) , 



1 



XI V W - Wi a/ {Wi - W'iliWi - W2) 

X {{w'i- Wi)(w4 - W2) E(2;; W) - {w^ - Wi){wi - W2) ^{z; W) 

Wi - W2 



+ {W4 - Wi){w3 + W4 + Wi+W2 + 2)Il 



'-■z-W 



-2 {W4 - Wi){w3 + 1){W2 + 1) n 

W4 {wi - W2) 



W4 — W2 

{W4 + l){Wi - W2) 



_{Wi + l){Wi - W2) ' ^' ^ 



+2 — wijw^ W2 n 



Wi [Wi - W2) 



(G.22) 
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where E, F and 11 are elliptic functions 

z = ArcSin 



W — Wi)(w4 — W2) 



W 



{w - W4){Wi - W2) 
(W3 - W4)(wi - W2) 



(G.23) 



(G.24) 



{W3 - Wi){w4 - W2) 

and Wi = e^' are the roots of uj'^{x) as given in the previous Section. The above expression 
evaluated at w = W2 yields 

UJ 



^/{Wl - W3){W4 - W2) 
+ {W4 - Wi) {W3 + W4 + Wi + W2 + 2)11 



{(W3 - Wi){w4 - W3) E{W) + {W4 - Wi) {W2 - W4) K{W) 
Wi - W2 



-2{W3 + 1){W2 + 1)(W4 - Wi) n 

+2 W3 W2 (.W4 — wi) n 



w 

W4 — W2 
{W4 + l){Wi - W2) 



_{Wi + 1){W4 - W2] 
W4 {Wi - W2) ^ 



■W 



(G.25) 



_Wi {W4 - W2, 

in which we used E(7r/2; q) =E(g), F(7r/2; q) =K(g) and n(p; 7r/2; q) = Il{p, q) 



G.2.3 Cubic interpolation 

In Section [0.21 we use the analytic approximation for a (a;) in Eq. fIG.lQp which involves 
the cubic interpolation 

a = Cq + Cix + C2X^ + C3X^ , (G.26) 

for niXi < X < K2 X2. The above coefficients Cj, with z = 0, . . . , 3, can be easily expressed 
in terms of the parameters A, B and D for the Morse potential flG.lSP as 

^2 \/D — 1 j K2^ X'l' (^2 X2 — Hi Xi) — 2 ^/D In ^2 -s/A^ Ki Xi (3 K2 X2 — Hi Xi) 



Co = Ki Xi 



2 ^/D [k2 X2 — Ki Xif 
{k2 X2 — Hi Xi) K2 X2 (2 Ki Xi + K2 X2) — 2 ^/D Ki Xi (4 k,2 X2 — k,i Xi) 



2 y/D {K2X2 - Ki Xi)^ 
12 y/D In ^2 VA^ Ki K2 Xi X2 



Co 



C. 



2 V -D {K2 X2 — KiXi) 

6 \/D In ^2 (/ti xi + K2 X2) + {2 \fD — 1 j {kx x\ + kx K2 x\X2 — 2 k-i X2 



(G.27) 



2 V D\Kx X\ — K2 X2) 
4 \/D In (^2 - (2 /D - 1) {k2 X2 - Kx xx) 



2 V-D (/€2 X2 - Kx Xx) 
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Finally, one can express everything in terms of the original parameters uo and 
use of Eqs. flG.15ap - (IG.15cl) . 



by making 
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